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SUMMARY 


Both  the  kinetics  and  the  mechanism  of  the  thermal 

oxidation  of  silicon  in  "dry"  oxygen  have  been  investigated 
18  15 

utilizing  the  0(p,a)  N nuclear  reaction.  The  oxygen-18 
concentration  profiles  were  determined  by  a trial  and  error 
deconvolution  procedure  of  the  recorded  alpha  spectra. 

The  kinetic  measurements  indicate  that  in  the  tem- 
perature range  investigated,  1Q89-5*K  to  1473-5#K,  the  para- 
bolic rate  constant  could  be  expressed  as 

kp  - 1(7. 5-3.0)  x 10"10  exp  - (~ ° 8^° ' Q9-*V)  ] 


It  was  found  that  there  was  good  agreement  between  the  silica 
scale  thickness  as  measured  by  the  single  spectrum  proton 
activation  technique  ajd  the  estimated  scale  thickness 
using  interference  colors. 

The  rate  controlling  process  in  the  thermal  oxidation 
of  silicon  in  the  temperature  range  1171-5‘K  to  1473-5#K  was 
determined  by  successive  isotopic  oxidation  studies  to  be  the 
transport  of  molecular  oxygen  through  the  adherent  silica 
scale.  The  exchange  process  between  the  permeating  oxygen  and 
the  network  oxygen  previously  reported  by  Ligenza  and  Spltzer 
was  confirmed.  It  was  found  that  the  kinetics  of  the  exchange 
and  permeation  processes  have  different  temperature  dependences 
and  that  the  kinetics  of  the  two  processes  are  approximately 
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equal  at  132S*K.  At  teaperatures  In  excess  of  132S°K  the 
exchange  kinetics  exceed  the  permeation  kinetics. 

The  combined  results  of  the  kinetic  and  mechanistic  studies 
confirm  the  theory  set  forth  by  Meek  to  explain  the  diffusion 
of  oxygen  in  vitreous  silica.  It  was  found  that  the  diffusivity 
governing  the  transport  of  molecular  oxygen  through  thermally 
formed  silica,  as  calculated  from  the  parabolic  rate  constant, 
was  in  agreement  with  the  diffusivity  measured  by  Morton  in  his 
permeation  experiment.  It  was  also  possible  to  calculate  an 
approximate  diffusivity  based  on  the  decrease  in  oxygen-18 
content  of  the  silica  scales  during  the  isotopic  oxidation 
experiments.  The  calculated  isotopic  diffusivity  was  found  to 
be  approximately  equal  to  the  diffusivity  measured  by  Williams 
In  his  exchange  experiments. 
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SECTION  I 


INTRODUCTION 

The  mechanism  of  the  thermal  oxidation  of 
silicon  Is  of  general  scientific  Interest  and  of 
specific  Interest  to  both  the  planar  device  Industry 
and  the  glass  Industry.  The  formation  of  a thin  layer 
of  silicon  dioxide  on  silicon  wafers  Is  of  importance 
In  the  masking  and  surface  passivation  of  planar  sili- 
con devices.  These  oxide  layers  can  be  formed  by  var- 
ious methods,  for  example  deposition  by  a vapor-phase 
reaction,  electro-chemical  oxidation  (anodization), 
plasma  reaction,  or  thermal  oxidation  by  one  of  the 
following  reactions: 

SI  + 02£S102 

Si  + 2H20^Si02  + 2H2 

In  Industrial  applications  silicon  dioxide  layers  are 
most  frequently  formed  by  thermal  oxidation. 

The  glass  industry  is  Interested  in  the  mechan- 
ism of  the  thermal  oxidation  of  silicon  In  the  hope 
that  It  will  lead  to  an  understanding  of  the  diffusion 
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of  oxygen  sed  silica.  An  understanding  of  oxygen 

diffusion  in  fused  silica  can  lead  to  a better  under- 
standing of  many  physical  properties  such  as  optical 
absorption,  viscosity,  and  electrical  conductivity.  In 
this  study  It  will  be  shown  that  the  thermal  oxidation 
of  silicon  Is  governed  by  the  diffusion  of  molecular 
oxygen  through  the  silica  scale. 

Because  the  adherent  oxide  scales  formed  on 
silicon  are  very  thin,  traditional  tools  for  studlng 
oxidation  processes,  such  as  weight  gain,  Inert  markers, 
and  optical  microscopy,  are  difficult  to  apply.  Elllp- 
sometry  has  proven  to  be  effective  In  making  kinetic 
& measurements,  but  Is  not  effective  In  determining  reac- 

tion mechanisms.  A modification  of  the  single  spectrum 

I proton  activation  technique  (1),  originally  developed  at 

Case  Western  Reserve  University  to  study  oxygen  diffu- 
sion in  oxide  ceramics,  was  used.  This  tecnnlque, 
because  of  Its  sub-micron  resolution,  allows  determi- 
nation of  both  the  kinetics  and  mechanism  of  the  reac- 
tion. 
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SECTION  II 


LITERATURE  SURVEY 

2 . 1 Theories  of  Oxidation 

2.1.1  Linear  Oxidation  Kinetics 

Linear  oxidation  kinetics  hold  when  the  rate 
of  oxidation  Is  constant,  or 


y * k£t  + constant 

where  k£  Is  the  linear  rate  constant,  y is  the  scale 
thickness  and  t is  time.  See  Figure  1.  An  example 
(2)  of  systems  exhibiting  this  type  of  kinetics  is  when 
the  oxidizing  environment  can  reach  the  metal  surface 
through  cracks  or  pores  in  the  oxide  coating. 

2.1.2  Parabolic  Oxidation  Kinetics 

If  the  oxide  scale  formed  is  protective,  a 
diffusion  process  through  the  scale  will  control  the 
oxidation  kinetics  and  parabolic  behavior  will  be 
observed.  Here  the  rate  of  scale  formation  will  be 
inversely  proportional  to  the  scale  thickness. 


y a 2kpt  + constant 

See  Figure  2. 

2.1.3  Linear-Parabolic  Oxidation  Kinetics 

The  1 Inear-parabol 1c  oxidation  kinetics  are 
derived  here  for  the  case  where  the  oxidation  reac- 
tion takes  place  at  the  metal-oxide  interface  (3,4) . A 
similar  derivation  could  be  presented  for  the  case 
where  metal  diffuses  through  the  scale  and  the  reac- 
tion takes  place  at  the  oxide-gas  Interface. 

If  oxidation  Is  to  occur,  the  Incoming  oxygen 
must  undergo  three  processes,  which  are  Illustrated  In 
Figure  3.  These  steps  are: 

1)  The  transport  of  the  oxygen  from  the  bulk  gas  to 
the  oxide-gas  Interface. 

2)  The  diffusion  of  the  oxygen  as  an  Ion,  atom,  or 
molecule  across  the  existing  oxide  layer. 

3)  The  reaction  of  the  oxygen  with  the  metal  at  the 
metal-oxide  Interface. 

The  fluxes  Jj,  J2,  and  J3  correspond  to  the 
three  processes  In  the  oxidation  process  and  can  be 
expressed  as  follows.  The  flux  of  oxygen  from  the  bulk 
gas  to  the  oxide  gas  Interface,  Jj,  Is  given  by: 
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J1  * K<ce  - c0> 


(3) 


The  proportionality  coefficient,  K,  is  the  surface- 
exchange  coefficient  and  Cg  and  CQ  are  given  In  Figure  3. 

The  flux  of  oxygen  through  the  existing  oxide 
layer,  J2,  Is  given  by  Fick's  First  Law: 

J2  = -8“  (4) 

ay 

where  s is  the  mobility.  If  the  solution  is  Ideal  then 
u **  RTln  C and  Equation  4 can  be  rewritten  as: 


aC 


If  the  Initial  oxide  scale  Is  thin  then  Fick's  First 
Law  can  be  approximated  as 

J2  ‘ (6) 

The  flux  corresponding  to  the  reaction  at  the  metal- 
oxide  Interface,  J3,  is  equal  to  the  eactlon  rate 
constant,  k$,  times  the  concentration  of  oxygen  at  the 
Interface: 


* 


By  definition  of  steady  state,  Jx  - J2  ■ J3  - J , 

It  Is  possible  to  solve  for  the  two  Interfacial  concen* 
tratlons  of  oxygen: 


1 + x + 1 r 


, k.y 
(1  + -$-)  c. 

D e 


1 + — =-  + 
K 


k5  x V 


As  D becomes  small  the  limiting  forms  of  Equa- 
tions 8 and  9 are  C,-*-  0 and  C^-*-  C and  the  oxidation  Is 

1 o e 

diffusion  controlled.  In  the  other  limit  as  0 becomes 
large  * CQ  * Ce/(1  + / K ) and  the  oxidation  process 

Is  reaction  controlled.  These  two  limits  are  shown 
schematically  In  Figure  4. 

dy 

The  rate  of  oxide  growth  Is  given  by  Nx — , where 

dt 

Ni  Is  the  number  of  oxygen  molecules  Incorporated  in  a 
unit  volume  of  oxide. 

At  steady  state 


Ni — ■ J » 


"sC. 

ks  k y 
1 + -=•  + — 

K D 


Applying  the  boundary  condition  that  y (0)  ■ y^, 
where  y1  Is  the  Initial  oxide  thickness,  and  solving 
the  differential  equation  the  following  Is  obtained 


where 


y ♦ Ay  « B(t  + t) 


1 ; 

A - 2D  ( — ♦ — ) 
k K 


(11) 


(12a) 


(12b) 


y,  ♦ Ay, 

T - ( ) 


(12c) 


Solving  for  the  oxide  thickness  as  a function 


of  time  yields 


A t ♦ t . A 

y ■ — <1  ♦ ( )>*  - - 

2 A2/4B  2 


Equation  11  reduces  to  linear  oxidation  kinetics 

2 

at  very  short  times,  l.e.,  (t  + i)  <<  A /4B.  The  linear 
rate  constant  kp  Is  given  by 


ksKCe 


* A (k  ♦ K)N, 


2 

For  long  times,  i.e.,  t>>A  / 4 B , Equation  11  reduces  to 
parabolic  oxidation  kinetics  with  the  parabolic  rate 
constant  being  given  by 


k = B 
P 

Figure  5 shows 
oxidation  kinetics. 

2.1.4  Wagner's  Theory  of  High  Temperature  Oxidation 

Carl  Wagner  developed  a theory  of  high  tempera- 
ture oxidation  (5-9)  In  which  the  oxidation  process  is 
diffusion  controlled,  therefore  leading  to  parabolic 
kinetics.  Wagner  assumes  that  positively  charged  metal 
Ions,  negatively  charged  nonmetal  Ions,  and  electrons 
are  the  only  mobile  species  In  the  oxide  and  that  the 

4 

diffusion  of  neutral  species  can  be  neglected.  If  this 
Is  the  case  there  are  two  limiting  cases:  (1)  metal 
Ions  and  electrons  diffuse  from  the  metal-oxide  inter- 
face to  the  oxide-gas  Interface  and  (2)  ele:trons  dif- 
fuse from  the  metal-oxide  Interface  to  the  oxide-gas 
Interface  while  nonmetal  Ions  diffuse  In  the  opposite 
direction.  If  the  diffusion  of  each  the  various  species 
Is  independent  of  the  others  their  drift  velocity  can 
be  expressed  as 


20C 


(15) 


a plot  of  the  linear-parabolic 
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its  two  1 


*1 ( 


id  x.edi 

J-  ♦ -1 ) 


N0dy 


(16) 


where  8 Is  the  mobility,  y the  chemical  potential,  N0 
Avogadro's  number,  z the  valence,  e the  electronic 
charge,  and  4 the  local  electrical  potential. 

The  number  of  moles  migrating  per  unit  cross 
section  per  unit  time  Is  then  given  as 


" c1v1  “ 


dy,.  z.ed* 

Ci  ( — ) (17) 

1 1 N„dy  dy 


where  c Is  the  concentration. 

Using  electrical  data  for  the  rates  of  Ion 
migration  In  terms  of  equivalents  per  unit  area  per 
second  Wagner  derived  the  following: 


n 300  y •'  1 1 

-f3-  - [( )N0e/  y ( -Htj  + t2)t3wdy  ]( ) (18) 

a 96500  y * | z 2 1 y aA y 


Here  y Is  the  chemical  potential  of  the  nonmetal,  y ' , 
the  chemical  potential  of  the  nonmetal  at  the  metal- 
oxide  Interface,  yy",  the  chemical  potential  of  the 
nonmetal  at  the  oxide-gas  Interface,  t the  trans- 
ference number,  k the  electrical  conductivity,  a the 
area,  and  the  subscripts  1,  2,  and  3 refer  to  the  metal 
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Ions,  nonmetal  Ions  and  electrons  respectively.  The 
part  of  Equation  18  In  square  braces  can  be  set  equal 
to  the  rational  rate  constant,  k.  for  the  oxidation 

r 

reaction. 

Assuming  the  transference  number  for  electrons, 
t^.  Is  approximately  unity  and  using  the  definition  of 
a transference  number  and  Equation  17,  one  can  solve  for 
the  electrlal  potential  gradient: 

d*  du3 

— - (19) 

dy  N0edy 

By  combining  Equation  17  and  19  It  Is  now  possible  to 
solve  for  the  number  of  moles  of  Ions  migrating  per 
unit  area  per  unit  time: 


hi  3iCx  dux  zxdu3 

— * ( ) ( + ) 

a Nq  dy  dy 


^2  dll2  Z2du3 

» ( ) (—  ♦ J—1) 

a N0  dy  dy 


(20a) 


(20b) 


If  chemical  equilibrium  Is  to  be  maintained  in 


the  system  then 


uM  * ui  + z j v 3 


(21a) 


>*n  = W2  + Z2V3 


(21b) 


where  the  subscripts  M and  0 refer  to  neutral  metal 
and  oxygen  atoms  respectively.  By  definition  of  the 


chemical  potential  In  terms  of  activities,  a, 


duM  ■ RTd(ln  aM) 


(22a) 


dyQ  = RTd(ln  aQ) 


(22b) 


Combining  Einstein's  Equation  relating  mobility 
and  the  dlffuslvlty  of  a species  with  Equations  20  and 
21  the  following  are  obtained: 

fix  din  a 


— - -DjCj ( 

a dy 


‘M 


) 


(23a) 


n2 


din  an 

-i , .d2C2( a 

a dy 


) 


(23b) 


If  the  oxide  Is  stolchometrlc  the  sum  of  the  dlf- 
fuslvltles  will  be  essentially  the  same  for  all  y's. 

If  Equations  23a  and  23b  are  multiplied  by  dy,  added, 
and  Integrated  over  the  thickness  of  the  oxide  the  fol- 
lowing Is  obtained: 


;2ini  z2n2 

( + )Ay 


= " 

SM 


-Zici/  Dxdln  aM  + |z2|c2/  D2dln  a. 
aM ' a0* 


(24) 
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The  activities  can  be  related  by 


| z 2 1 d 1 n aM  ♦ Z}  din  aQ  ■ 0 
and  the  rational  rate  constant  expressed  as  either 


V zi0x 

kr  “ C.n/  ( 

r *q  *0'  I'll 


♦ 0^)  din  a( 


(26a) 


M I z 2 I D2 

kr  ' .(D*  * 1 d,n  *M 


(26b) 


where  c#q  » c ^ - c2|z2|. 


If  the  oxide  Is  nonstolchometrlc  deviations 
can  occur  In  four  different  ways.  Excess  metal  can  be 
accounted  for  by  either  cation  Interstitials  and  free 
electrons  or  anion  vacancies  and  fr».*e  electrons.  A 
metal  deficit  can  be  accounted  for  by  either  cation 
vacancies  and  electron  holes  or  anion  Interstitials 
and  electron  holes. 

For  the  case  of  excess  metal  Ions  the  rational 
rate  constant  can  be  expressed  as 


kr  • <»  ♦ 'll  C,,"!’  (1  - 


CM1)' 


where  the  subscript  (1)  refers  to  the  Interstitial 


c ^ 

position  and  kT. 


If  there  is  a metal  deficit,  the  rational  rate 
constant  Is 


p02‘  l/8 

kr  * (1  + zi)ceqDi"n  - ( ) ] (28) 


p02" 


£ 1 H 

where  Oj " ■ ( LLljgj  kj  and  the  subscript  (v) 

Cj  v ' 

refers  to  vacant  sites. 

For  the  case  of  excess  oxygen  Ions  the  rational 
rate  constant  Is 


p0  • / 

kr  3 (1  + |22l)ceq02'[l  - (~^-)  3]  (29) 


p0," 


c2 " 

where  02  • « ( ^)32(ijlc7. 


In  the  case  of  an  oxygen  deficit  the  rational 
rate  constant  can  be  expressed  as 


P0,\  Vs 


kr  " (1  + ) z2 I )ceq°2 "[1  - ( ) ] (30) 


p0, 


2.1.5  logarithmic  Oxidation  Kinetics 

The  logarithmic  rate  law  can  be  stated  as  follows 

1 

- * Y - Z log  t (31) 

y 

where  Y and  Z are  related  to  the  oxidizing  pressure* 
the  oxidizing  temperature,  a characteristic  distance 
dependent  on  the  contact  potential  difference  (work 
function),  and  temperature,  and  a limiting  thickness. 
See  Figure  6. 

Oxidation  kinetics  of  this  type  are  observed  for 
many  metals  at  low  temperatures  (e.g.  aluminum  at  room 
temperature).  A possible  mechanism  for  this  type  of 
kinetics  has  been  set  forth  by  Cabrera  and  Mott  (10- 
12).  It  Is  assumed  that  oxygen  atoms  are  adsorbed  on 
the  surface  of  the  oxide,  which  produces  vacant  energy 
levels  for  electrons.  Electrons  can  nc v pass  through 
the  oxide  layer  from  the  metal  to  the  oxygen  atoms  by 
either  thermionic  emission  or  tunnelling,  to  equalize 
the  surface  state  levels  and  the  Fermi  energy  of  the 
metal.  The  electrons  will  then  combine  with  the  oxy- 
gen atoms  at  the  surface,  setting  up  a potential  dif- 
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2.2  Silicon  Oxidation  K1 netlcs 

The  kinetics  of  the  thermal  oxidation  of  silicon 
In  dry  oxygen  obeys  a linear-parabolic  rate  law  (3,4, 
13-17)  with  the  Initial  linear  portion  corresponding 
to  a thickness  of  approximately  230  A (3,4). 

The  activation  energy  for  the  linear  portion  of 
the  thermal  oxidation  reaction  is  1.98±0.02  eV  (3,4, 
17).  However,  there  is  a large  discrepancy  In  the 
reported  activation  energy  for  the  parabolic  portion 
of  the  reaction.  The  experimentally  determined  acti- 
vation energies  for  this  portion  of  the  reaction  fall 
In  two  ranges,  1.27±0.06  eV  (3,4,14-16)  and  1.80+0.10 
eV  (17-19).  The  writer  can  see  no  correlation  between 
the  activation  energy  determined  and  any  of  the  experi- 
mental parameters  (e.g.  temperature  range,  pressure 
range,  crystal  orientation,  or  type  of  oxidation). 

It  has  been  observed  by  Irene  (17)  that  both  the 
linear  and  parabolic  rate  constants  Increase  with 
increasing  water  content  In  the  ambient  gas.  See  Table 
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TABLE  1 

THE  LINEAR  AND  PARABOLIC  RATE  CONSTANT 
EXPRESSED  AS  A FUNCTION  OF  WATER 
IN  THE  AMBIENT  AND  ORIENTATION  ACCORDING 
TO  IRENE  (17^ 


TEMPER- 

ATURE 

(°C) 

ORIEN 
TAT ION 

h2o 

in  02 
(PPM) 

la 

(A/min) 

STAN 

CARD 

DEVI 

AT  ION 

STAN 

DARD 

2 kp  DEVI 

(A  / min)  ATION 

800 

(100) 

<1 

0.67 

0.04 

310 

30 

*25 

0.83 

0.02 

430 

20 

(no) 

<1 

1.3 

0.10 

480 

50 

*25 

1.3 

0.03 

800 

40 

(111) 

<1 

1.4 

0.05 

690 

20 

*25 

1.3 

0.03 

1,200 

60 

927 

(100) 

<1 

5.3 

0.20 

2,700 

100 

*25 

5.3 

0.30 

4,300 

260 

(no) 

<1 

11 

0.90 

3,600 

150 

*25 

8.3 

0.60 

5,300 

250 

(111) 

<1 

8.9 

0.20 

5,300 

110 

*25 

9.7 

0.60 

6,900 

290 

996 

(100) 

<1 

16 

0.80 

7,700 

240 

~25 

13 

1.0 

12,000 

890 

(no) 

<1 

29 

3.4 

9,400 

360 

*25 

28 

0.80 

12,000 

210 

(111) 

<1 

26 

2.0 

11,000 

370 

*25 

28 

0.80 

13,000 

230 

22 


1 

It  has  been  reported  by  Deal  (15)  and  Irene  (17) 
that  there  Is  an  orientation  effect  on  the  thermal 
oxidation  of  silicon.  It  can  be  seen  from  Irene's 
(17)  data  In  Table  1 that  the  dependence  of  the  linear 
rate  constant  on  orientation  Is 

(110)  > (111)  > (100) 

Irene's  rationalization  of  this  Is  that  the  order  and 
the  magnitude  of  the  rate  constants  are  In  agreement 
with  the  number  of  silicon  to  silicon  bonds  In  each 
of  these  planes. 

Irene  (17)  also  reports  that  there  is  an  orienta- 
tion dependence  of  the  parabolic  rate  constant  which  Is 

f> 

(111)  > (110)  > (100) 

To  date  this  dependence  has  not  been  explained  since 
the  orientation  of  the  silicon  should  not  effect  the 
diffusion  of  oxygen  through  the  silica  scale  forming  on 
the  surface  of  the  crystal. 

It  has  been  reported  by  Oeal  (15)  that  the  kine- 
tics of  thermal  oxidation  of  silicon  Is  Independent  of 
Impurity  type  (phosphorus  or  boron  dopants)  and  also 
of  impurity  concentration  In  the  range  of  0.001  to  50 
ohm-centimeter.  He  also  found  that  the  properties 
of  the  oxide  formed  were  Independent  of  both  Impurity 
type  and  concentration. 
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2.3  S11 Icon  Oxidation  Mechanisms 


Because  of  the  parabolic  behavior  of  the  thermal 
oxidation  of  silicon  It  Is  generally  agreed  that  the 
reaction,  at  least  at  long  times  and  high  temperatures , 

Is  controlled  by  the  diffusion  of  either  silicon  Ions, 
oxygen  Ions,  oxygen  atoms,  or  oxygen  molecules  (electron 
diffusion  Is  not  considered  here  because  of  the  corre- 
sponding high  transport  number)  through  the  silica  scale. 
If  we  look  at  the  activation  energy  of  the  parabolic 
portion  of  the  rate  constant,  we  find  agreement  between 
the  activation  energy  of  1.27  electron  volts  (3,4,14-16) 
and  the  activation  energy  for  oxygen  diffusion  In  fused 
silica  determined  by  Williams  (20)  and  Norton  (21)  (See 
Section  2.4).  This  agreement  Indicates  the  possibility 
that  the  thermal  oxidation  of  silicon  is  controlled  by 
oxygen  diffusion  assuming  these  are  the  correct  acti- 
vation energies  for  oxidation  and  diffusion. 

Studies  of  the  oxygen  partial  pressure  dependence 
of  the  thermal  oxidation  reaction  have  shown  the  linear 
rate  constant  to  be  pressure  Independent.  It  has  been 
reported  by  Atalla  (18)  that  the  parabolic  rate  con- 
stant varies  as  the  partial  pressure  of  oxygen  to  the 
four-fifths  power;  however,  Llgenza  and  Spitzer  (22) 
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and  Deal  and  Grove  (4)  reported  the  parabolic  rate  con- 
stant to  be  directly  proportional  to  the  oxygen  partial 
pressure.  This  dependence  Is  the  same  as  that  found 
for  oxygen  diffusion  In  fused  silica  (21)  and  Indicates 
that  transport  occurs  by  the  diffusion  of  neutral  oxygen 
molecules. 

Further  confirmation  that  thermal  oxidation  occurs 
by  the  transport  of  oxygen  rather  than  silicon  was  pre- 
sented by  Atalla  (18).  Silicon  samples  were  partially 
oxidized  and  then  a layer  of  radioactive  antimony 
evaporated  on  the  thermal  silica  scale.  The  oxidation 
was  then  continued  and  the  samples  analyzed  for  anti- 
mony. The  radioactive  antimony  was  found  at  the  oxide- 
gas  interface  (See  Figure  7)  Indicating  oxygen  diffusion 
to  be  the  rate  controlling  process  In  the  thermal  oxida- 
tion of  silicon. 

Further  confirmation  that  the  oxidation  reaction 
occurs  at  the  silicon-silica  Interface  was  provided  by 
the  experiments  of  LI  jenza  and  Spitzer  (22).  In  their 
experiments  silicon  coupons  were  oxidized  successively 
In  oxygen-16  and  oxygen-18  and  the  Isotopic  shift 
determined  by  the  fundamental  Infrared  absorption  bands 
of  silica  which  show  appreciable  Isotopic  shift.  See 
Figure  8. 


(A)  Initial 


♦ 950A-#- 


Figure  7:  The  Initial  and  final  distributions  of  radio- 
antimony found  by  Atatta  (18)  at  1193°K. 
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Figure  8:  The  fundamental  infrared  absorption  bands 
of  silica  which  show  appreciable  isotopic  shift.  Curve  (1) 
is, for  a sample  of  S1180?t  and  curve  (2)  is  for  a sample  of 
Si 1602. 
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Llgenza  and  Spltzer  found  that  If  they  oxidized  a 
silicon  coupon  In  oxygen-16  and  then  continued  the 
oxidation  In  oxygen-18,  almost  the  entire  oxide  scale 
was  composed  of  oxygen-18.  If  the  Isotopic  order  was 
reversed,  they  found  the  silica  to  be  almost  entirely 
oxygen-16.  It  becomes  apparent  from  these  results 
that  not  only  does  growth  occur  during  the  second  oxida- 
tion, but  also  exchange  of  the  ambient  oxygen  with  the 
oxygen  of  the  silica  scale.  The  species  responsible 
for  the  oxidation  and  exchange  reactions  must  be  oxygen 
since  such  exchange  cannot  be  accounted  for  by  the  dif- 
fusion of  silicon  or  vacancies. 

It  has  been  found  by  Jorgensen  (13)  that  the 
kinetics  of  the  thermal  oxidation  reaction  of  silicon 
was  controlled  by  the  migration  of  an  oxygen  species 
and  the  kinetics  could  be  effected  by  the  application 
of  an  electric  field.  Jorgensen  first  performed  a mod- 
ified marker  experiment  In  which  he  sputtered  a thin 
strip  of  platinum  onto  a silicon  wafer.  A layer  of 
silicon  monoxide  was  evaporated  on  top  of  the  silicon 
and  platinum  and  then  oxidized  to  silicon  dioxide. 

(See  Figure  9 for  schematic  of  sample)  The  entire 
assemblage  was  then  oxidized  and  the  platinum  marker 
located.  It  was  found  that  the  platinum  was  not  located 


at  the  silicon-silica  Interface,  Implying  that  oxidation 

occurs  by  the  migration  of  oxygen. 


Pt 


. 


Figure  9.  A schematic  representation  of 

Jorgensen's  (13)  modified  marker 
experiment. 


In  a second  experiment  Jorgensen  concluded  that  the  rate 
controlling  species  in  the  oxidation  reaction  was  oxygen  ions. 
If  the  kinetics  are  controlled  by  the  diffusion  of  a charged 
species  then  it  should  be  possible  to  influence  the  reaction 
rate  by  the  application  of  an  electric  field,  while  if  a neu- 
tral species  were  migrating  an  electric  field  should  have  no 
effect  on  the  kinetics.  Jorgensen  found  that  if  he  made  the 
oxide-gas  interface  negative  with  respect  to  the  metal oxide 


interface  the  reaction  rate  was  increased,  while  if  he  made 
the  oxide-gas  interface  positive  with  respect  to  the  metal- 
oxide  interface  he  could  either  retard  or  completely  stop 
the  oxidation. 


Jorgensen's  experiment  has  been  reinterpreted 
by  Raleigh  (23)  In  terms  of  the  electrochemistry  of  the 
process.  A schematic  of  the  apparatus  used  In  Jorgen- 
sen's experlement  Is  shown  In  Figure  10.  A porous 
platinum  electrode  was  attached  to  each  end  of  a par- 
tially oxidized  single  crystal  of  silicon  and  the 
assembly  heated  In  an  oxygen  atmosphere.  The  rate  of 
silica  growth  at  each  end  of  the  sample  was  the  moni- 
tored in  the  presence  of  an  electric  field  (approxi- 

4 

mately  10  volts  per  centimeter  across  the  silica)  and 
in  the  absence  of  an  electric  field.  Raleigh  argues 
that  the  external  application  of  an  electric  field  by  itself  can- 
not provide  a steady-state  driving  force  for  the  oxidation  reaction. 
This  is  because  in  any  diffusion-controlled  oxidation,  the  net 
process  is  required  to  be  the  transport  of  an  electroneutral  species 

(this  does  not  rule  out  the  transport  of  ions  as  long  as  charge 
neutrality  is  maintained)  through  the  growing  oxide  scale. 


i 


He  also  argues  that  oxygen  diffusion  is  the  result  of  a chemical 
potential  gradient  which  is  fixed  on  both  sides  of  the  film  and 
cannot  be  effected  by  an  externally  applied  electric  field. 

Raleigh  explains  his  argument  by  assuming  that  oxygen 
ions  are  transported  across  the  growing  silica  film.  In  order 
to  maintain  electrical  neutrality  the  transport  of  oxygen  ions 
must  be  compensated  for  by  the  transport  of  an  equal  amount  of 
electronic  charge  in  the  opposite  direction.  Therefore  the 
field  used  to  accelerate  the  ionic  transport  will  retard  the 
electronic  transport  and  the  result  of  this  unbalanced  trans- 
port will  be  the  transport  of  unbalanced  charge.  In  the 
steady-state  the  uncompensated  charge  will  set  up  a field  equal 
and  opposite  to  the  applied  field  which  will  tend  to  cancel  it 
and  return  the  system  to  the  field  state. 
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Raleigh  concludes  that  the  reason  Jorgensen 
could  Influence  the  oxidation  rate  by  the  application 
of  an  electric  field  was  in  the  use  of  a battery  to 
apply  the  electric  field.  The  battery  circuit  pro- 
vided an  external  circuit  for  electronic  flow,  oxygen 
ion  transport  then  required  only  Ionic  conductivity, 
allowing  what  may  have  been  only  a minor  transport 
mechanism  because  of  a low  electronic  conductivity 
to  assume  a major  significance.  This  in  essence 
set  up  an  electrolysis  cell. 


2.4  Oxygen  Diffusion  in  Fused  Silica 

There  is  a large  variation  in  the  di ffusi vi ti es 
of  oxygen  reported  in  the  literature  (20,  21,  24,  25). 
The  reported  values  of  the  pre-exponential  term  vary 
over  seven  orders  of  magnitude  and  the  reported  activa- 
tion energies  range  from  1.17  electron  volts  per  atom 
to  3.07  electron  volts  per  atom.  See  Table  2 and  Figure 
11. 

The  two  most  commonly  quoted  values  for  oxygen 
diffusion  in  silica  are  those  given  by  Williams  (20): 


9 1.26 

D = 2.0x10"  exp( ) (32) 

kT 

and  Norton  (21): 

4 1.17 

D = 2.88x10"  exp( ) (33) 

kT 
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TABLE  2.  — Oxygen- 18  dlffuslvitles  in  fused  silica. 


Investigator 


Temperature 

Range 

(°K) 


— X 

(cm2/?ec)  (kcal/mole) 


1.  Haul,  Dumbgen 

1173-1523 

4.3  x 10'6 

56 

2.  Norton 

1173-1373 

2.88  x 10‘4 

27 

3.  Sucov 

1198-1498 

1.5  x 10'2 

71 

4.  Williams 

1123-1523 

2.0  x 10‘9 

29 

Meek  (26)  has  proposed  a model  that  can  bring  these 
two  experimental  values  Into  agreement. 

Meek  first  looks  at  the  experiment  preformed 
by  Williams  and  Norton.  Williams'  (20)  experiment 
Involved  measuring  the  decrease  of  oxygen-18  con- 
centration In  the  ambient  gas  In  contact  with  either 
silica  fibers  or  bulbs.  He  found  that  the  dlffuslvlty 
was  proportional  to  the  partial  pressure  of  oxygen. 

1 
j j 

ular  transport. 

Norton  (21)  on  the  other  hand  determined  the 

I 

dlffuslvlty  of  oxygen  through  silica  membranes  by 
measuring  the  permeation  of  oxygen-16  In  silica.  Here 
he  applied  a positive  oxygen  pressure  to  one  side  of 
a membrane  and  followed  the  mass  32  peak  with  a mass 
spectrometer  connected  directly  to  the  other  side  of 
the  membrane. 

The  activation  energies  reported  by  both  Williams 
and  Norton  agree  within  experimental  error  and  are 
nearly  the  same  as  the  activation  energy  (1.15  eV)  for 
diffusion  or  argon  In  fused  silica,  wnich  has  the  same 
size  an  oxygen  molecule  (27).  The  result  leads  on  to 
believe  that  the  diffusing  species  are  neutral  oxygen 

molecules  which  pass  through  Interstlcies  In  the 
silica  network. 


36 


— 


Meek  feels  that  the  tracer  dlffuslvity  as  measured 
by  an  exchange  experiment,  such  as  Williams  performed. 

Is  actually  the  dlffuslvity  times  the  fraction  of  time 
the  tracer  Is  Incorporated  In  the  diffusing  species, 
providing  the  exchange  reaction  Is  rapid  and  diffusion 
Is  the  rate  limiting  process.  The  dlffuslvity  as  measured 
by  Williams  Is  then  an  effective  dlffuslvity,  which 

Is  given  by 

I 

°ef f * 0 ((?'■;  -c' ) <34) 


or 


D 


+ C' 

) 

C* 


(35) 


I 

where  C Is  the  Interstitial  solubility  of  oxygen  In 

N 

silica  and  C Is  the  concentration  of  oxygen  Incorporated 
In  the  silica.  Equations  34  and  35  assume  that  diffusion 
Is  the  rate  limiting  step  and  exchange  between  the  diffus- 
ing oxygen  and  the  network  oxygen  Is  rapid.  This  exchange 
reaction  can  also  be  used  to  explain  the  results  previ- 
ously reported  by  Hutchison  (28)  and  Ligenza  and  Spitzer 
(22). 
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2.5  Charged  Particle  Activation  Analysis  by  D1 rect 
Nuclear  Reactl ons 
2.5.1  General 

The  Identification  of  light  elements  (Z<^8), 
which  have  been  Introduced  In  samples  as  tracers. 
Impurities,  and  dopants.  Is  easily  accomplished  by 
Inducing  nuclear  reactions  by  charged  particles 
(protons  and  deuterons).  This  type  of  analysis  has 
been  finding  increasing  use  In  the  past  ten  years  to 
study  diffusion  (1,  29-34),  oxidation  (35),  electro- 
chemical phenomena  (36-39),  etc.  These  procedures 
have  the  following  properties  (40): 

(a)  Nuclear  microanalysis  Is  restricted  to  the 
first  few  microns  of  the  sample. 

(b)  With  most  of  the  reactions  there  Is  no 
natural  background. 

(c)  Most  of  the  reactions  have  a high  Q value* 

(see  Table  3)  which  allows  the  use  of  low 
energy  particles,  thus  reducing  local  heat- 
ing of  the  target,  and  avoiding  competing 
nuclear  reaction  (See  (d)). 

(d)  Coulomb  barriers  prevent  the  reaction  of  most 
nuclei  larger  than  chlorine  with  the  low  energy  particles 

* • The  Q value  of  a reaction  Is  the  energy  equlv- 
of  the  mass  difference  between  the  reactant  and  pro- 
duct nuclei . 
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Isotope 

Qo 

(Hey) 

Isotope 

Qo 

(MeV) 

Isotope 

Qo 

(MeV) 

17.347 

(p.a)  reaction*: 

*L1  4.02 

9 

„Be 

2.125 

“r 

8.582 

“o 

3.970 

31 

P 

1.917 

19 

F 

8.119 

Cl 

3.030 

27 

Al 

1.594 

IS 

N 

4.S64 

2) 

Ns 

2.379 

17 

io° 

1.197 

10 

B 

17.819 

(d.a)  reaction*: 

' 8 8.022 

10B 

32 

S 

1.147 

4.890 

6u 

22.36 

*8 

7.693 

“o 

4.237 

7n 

14.163 

*Be 

7.152 

30S1 

3.121 

14n 

N(-o) 

13.579 

*5Mg 

7.047 

“o 

3.116 

**F 

10.038 

23n* 

6.909 

“Mg 

2.909 

^0 

9.812 

27ai 

6.701 

24Mg 

1.964 

"(«,) 

9.146 

”« 

6.012 

•V, 

1.421 

MP 

8.170 

nc 

5.167 

“c 

<0 

*°B 

9.237 

(d.p)  reaction*: 

* 0 5.842 

4.212 

“MS 

8.873 

5.499 

“c 

2.719 

%,) 

8.615 

*Vg 

5.106 

1#o 

1.919 

2,si 

8.390 

•u 

5.027 

,Bo 

1.731 

32s 

"si 

6.418 

23n* 

4.734 

nN(Ps) 

B 

1.305 

6.253 

*Be 

4.585 

1.138 

”c 

5.947 

19f 

4.379 

15h 

0.267 

JIP 

5.712 

30si 

4.367 

7li 

<0 
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employed.  This  allows  the  determination  of  light 
elements  In  heavy  substrates. 

(e)  The  nuclear  reactions  employed  are  very  spe- 
cific, In  that  two  Isotopes,  even  of  the  same 
element,  react  quite  differently. 

(f)  The  results  obtained  are  Independent  of  the 
matrix  In  which  the  observed  Isotope  Is  embed- 
ded, since  the  physical  or  chemical  state  of 
the  nuclei  does  not  effect  the  reaction  yield. 

(g)  The  results  are  quantitative. 

(h)  The  method  Is  essentially  nondestructive. 

(1)  The  sensitivity  Is  high  In  that  as  little  as 
-12 

10  grams  of  matter  per  square  centimeter 
can  be  determined  under  Ideal  conditions. 

(j)  Concentration  versus  depth  relationships  to 
tens  of  microns  can  be  determined  with  sub- 
micron resolution  by  one  of  three  methods: 

(1)  serial  sectioning  (30,  31,  34,  35),  (2) 
resonance  techniques  (32,  41-44),  and  (3) 
single  spectrum  method  (1,  30,  31,  35,  40). 

2.5.2  Single-spectrum  Proton  Activation  Analysis  Utlli- 
16  is 

zing  the  0(p,o)  N Reaction  (1,44) 

If  a solid  containing  a concentration  gradient  of 
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oxygen-18  Is  Irradiated  with  a monoenergetlc  beam  of 
protons  and  the  energetic  alpha  particles  simultane- 
ously detected,  It  Is  a stral ghtforward  process  to 

convert  the  alpha  spectrum  to  a concentration  versus 

18  15 

depth  profile.  Since  the  0(p,o)  N reaction  Is 
exoerglc  with  a Q value  of  3.98  MeV  (see  Table  3), 
each  energetic  alpha  particle  emitted  can  be  associated 
with  a given  depth  In  the  sample.  A proton  Incident 
upon  a sample  with  energy  Ep(0)  will  lose  energy  almost 
linearly  as  It  traverses  the  sample  and  therefore  at  a 
depth  x from  the  sample  surface  will  have  an  energy 
Ep(x).  At  some  depth,  x,  the  proton  will  react  with 
an  oxygen-18  nucleus  producing  an  energetic  alpha  par- 
ticle of  energy  Ea(x,x)  which  depends  on  Ep(x)  in  the 
following  way  (45): 


(M  M E ( x ) ) **  M M E (x)  2 

E„(x,x)  3 [ — B— ° P Cos  i ± — 2 Cos  + 


Ma  + MN 


Ma  * V 


V + ep(x)(mn  " V „ 2 
}*] 

Ma  + 


(36) 


Here  Mp  Is  the  mass  of  a proton,  Ma  the  mass  of  an  alpha 
particle,  M^  the  mass  of  the  residual  nitrogen-15  nucleus, 
and  t|>  the  laboratory  angle  of  detection.  The  alpha  par- 
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tide  will  also  be  attenuated  as  It  travels  out  of 

the  sample  and  at  the  sample  surface  will  have  energy 

E (x,0).  As  can  be  seen  from  Figure  12,  there  will 

be  a one-to-one  correlation  between  alpha  energy  and 

depth,  If  one  knows  the  stopping  power  of  the  sample 

for  protons  and  alphas,  and  a concentration  versus 

depth  profile  Is  determined  from  the  number  of  alpha 

particles  of  each  energy  detected.  The  number  of 

2 

counts  In  a given  channel  , N.. , where  1 refers  to  the 
channel  number.  Is  given  by 

Nj  * Ltc(x.j  )o(xj  )axj an  (37) 

Here  L Is  the  volume  number  density  of  oxygen  atoms  In 
the  sample,  t the  Integrated  proton  flux,  C(x<)  the 
concentration  In  oxygen-18  fraction  at  mean  depth  x^, 
o(xi)  the  differential  cross  section  at  mean  depth  x^ 
(see  Figure  13  for  o(Ep)),  ax^  the  depth  Interval  with 
mean  depth  x^,  and  An  Is  the  solid  angle  of  observation. 

The  number  of  counts  determined  by  Equation  37 
holds  only  for  an  Ideal  spectrum  where  there  is  no  dis- 
tortion. In  actuality  the  spectrum  recorded  will  be 


2,A  channel  Is  a segment  of  the  multi-channel 
analyzer  which  records  and  stores  the  number  of  alpha 
particles  detected  In  a given  energy  Interval. 


mb/ste 


Ep  (KeV) 


Figure  13:  Absolute  cross  section  of  180(p,o)15N 
reaction  as  a function  of  incident  energy  (\j>=165°). 
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convoluted  with  a spreading  function,  P (Ej,Ej), 
which  Is  the  probability  that  an  alpha  particle 


In  channel  j.  The  distorted  spectrum,  [N^],  will 
then  be  given  by 


[Nj  ] 


5 Y<Ei-V 


(38) 


SECTION  III 
MECHANISTIC  MODELS 


3. 1 Introduction 

The  kinetics  of  the  thermal  oxidation  reaction 

Si  + 02  * Si02 

will  be  controlled  by  either  the  outward  diffusion  of 
silicon  ions  or  the  inward  diffusion  of  oxygen  atoms, 
molecules,  or  ions.  In  order  to  determine  the  reaction 
mechanism  it  is  necessary  to  perform  a marker  study. 

Conventional  marker  studies  of  the  Smiqelskas  and 
Kirkendall  (46)  type  are  difficult  to  perform  in  systems 
where  there  is  good  oxidation  resistance  due  to  the 
length  of  time  required  to  grow  an  oxide  scale  of  suf- 
ficient thickness.  It  has  been  estimated  by  Jorgensen 

s 

(13)  that  this  experiment  on  silicon  would  require  10 
years  using  a wire  marker  0.002  inch  in  diameter. 

Another  possible  type  of  marker  study  is  to  mark 
the  sample  surface  with  implanted  atoms  (47-48).  This 
method,  however,  can  lead  to  erroneous  results,  since  it 
must  be  assumed  that  such  "marker"  atoms  are  truly  inert. 
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1.e.,the  Implantation  process  does  not  affect  the  oxida- 
tion reaction,  and  the  Inert  gas  Is  Immobile. 


The  best  way  to  study  the  oxidation  of  an  oxida- 
tion resistant  material  Is  to  perform  a tracer  study 
similar  to  those  performed  by  Bardeen,  Brattaln,  and 
Shockley  (49)  and  Holt  and  Himmel  (50).  In  this  type  of 
study,  an  Isotope  of  one  of  the  reacting  species  is 
Incorporated  In  the  growing  scale.  By  locating  the 
Isotope  after  the  completion  of  oxidation  it  Is  possible 
to  determine  the  reaction  mechanism.  Holt  and  Himmel 
(50)  have  presented  the  expected  tracer  distributions  for 
the  case  where  oxidation  occurs  by  the  transport  of  an 
Ionic  species  and  either  the  cation  or  anion  Is  used  as 
a tracer.  These  profiles  are  shown  schematically  In 
Figure  14. 

Figures  14a  and  b show  the  expected  tracer  dis- 
tribution If  the  oxidation  reaction  is  controlled  by  the 
transport  of  metal  Ions  from  the  metal-oxide  Interface 
to  the  oxide-gas  Interface  with  the  reaction  occurring 
at  the  oxide-gas  Interface.  Figure  14a  shows  the  tracer 
distribution  If  a radioactive  cation  Is  used  as  the 
tracer  and  Figure  14b  the  tracer  distribution  If  a radio- 
active anion  Is  used  as  the  tracer.  Figures  14c  and  d 
show  the  corresponding  situation  for  the  case  where  the 
oxidation  reaction  is  controlled  by  the  transport  of  an 
Ionic  oxygen  species  and  the  chemical  reaction  occurs  at 


I 
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Cation  Mobile 


* . 14:  Expected  cation  (M*) 
and  anion  (X*)  tracer  distributions  in  the 
scale  under  conditions  (a)  and  (b),  res- 
pectively, where  the  cation  is  the  diffus- 
ing component.  The  corresponding  distribu- 
tions In  situations  where  the  anion  is  the 
mobile  species  are  Indicated  by  (c)  and 
(d),  respectively  (50). 
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the  metal-oxide  Interface.  Figure  14c  shows  the  case  If 
an  Isotope  of  the  metal  Is  used  as  the  tracer  and  Figure 
14d  the  case  If  an  oxygen  Isotope  Is  used  as  the  tracer. 

The  profiles  presented  in  Figures  14b  and  d cover 
the  situation  where  the  thermal  oxidation  of  silicon  Is 
controlled  by  the  transport  of  an  ionic  species  through 
the  adherent  silica  scale,  but  do  not  cover  the  possibil- 
ity that  the  thermal  oxidation  Is  controlled  by  the  tran- 
sport of  molecular  oxygen  through  the  silica  scale  as 
predicted  in  the  literature.  The  problem  Is  further  com- 
plicated If  the  oxidation  reaction  Involves  an  exchange 
reaction  along  with  the  transport  process  as  predicted  by 
Llgenza  and  Spitzer  (22). 

The  remainder  of  this  chapter  will  be  devoted  to 
a discussion  of  the  expected  oxygen-18  profiles  assuming 
the  oxidation  reaction  is  controlled  by  the  transport  of 
molecular  oxygen,  Ionic  oxygen,  and  ionic  silicon.  Each 
profile  has  been  developed  by  assuming  the  sample  has 
been  oxidized  In  pure  oxygen-18  to  produce  a scale  of  half 
the  desired  final  thickness.  The  ambient  gas  is  then 
changed  to  pure  oxygen-16  and  the  oxidation  continued  until 
the  final  thickness  is  reached.  The  sample  is  then 
analyzed  and  the  oxygen-18  concentration  versus  depth 
profile  obtained.  These  profiles  are  shown  as  oxygen-18 
concentration  as  a function  of  scale  thickness,  y,  with 
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y ■ 0 being  the  oxide-gas  Interface,  y » yi  the  Initial 
oxygen-18  scale  thickness,  and  y « 2y^  the  final  scale 
thickness. 

3.2  Oxidation  Controlled  by  the  Transport  of  Molecular 
Oxygen 

3.2.1  General 

If  the  thermal  oxidation  of  silicon  Is  controlled 
by  the  transport  of  molecular  oxygen  from  the  silica-gas 
Interface  to  the  silicon-silica  Interface  and  the  trans- 
port process  Is  similar  to  that  proposed  by  Meek  (26)  for 
oxygen  diffusion  In  silica  (the  permeating  oxygen  exchanges 
with  the  network  oxygen)  then  three  types  of  possible  pro- 
files exist  depending  on  the  relative  rates  of  permeation 
and  exchange.  In  the  first  case  the  rate  of  permeation 
Is  much  greater  than  the  rate  of  exchange  (the  case  for 
no  exchange  Is  a special  case  of  this);  the  second  case 
Is  for  the  permeation  rate  to  be  less  than  the  exchange 
rate;  and  the  third  case  is  Intermediate  to  the  first  two. 
Each  of  these  cases  will  be  discussed  below. 

Although  a mathematical  solution  for  these  three 
cases  will  not  be  presented,  it  should  be  possible  to 
express  the  permeation  rate  by  the  following  equation  (51): 

DSPa 

Permeation  Rate  * (39) 

y 

where  S is  the  solubility  of  oxygen  In  silica,  P the 
pressure  across  the  oxide,  and  a the  cross  sectional  area. 
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Two  exchange  reactions  will  also  occur.  The  first  ex- 
change occurs  at  the  sample  surface  and  involves  an  ex- 
change between  the  ambient  oxygen  and  the  surface.  This 
exchange  rate  can  be  expressed  by  the  following  equation: 

Surface  Exchange  Rate  ■ k (Ce  ” C0^  (40) 

where  K is  the  surface  exchange  coefficient,  Cfi  the  con- 
centration of  oxygen-18  in  the  ambient  and  CQ  the  concen- 
tration of  oxygen-18  at  y = 0.  It  has  been  assumed  in 
this  section  that  the  surface  exchange  coefficient,  K,  is 
equal  to  infinity.  The  other  exchange  occurs  between  the 
permeating  oxygen  and  the  network  oxygen;  this  exchange 
rate  can  be  expressed  by. 

Exchange  Rate  = K1([16023  - [Si1802])a'  (41) 

where  Kj  Is  the  exchange  coefficient,  [ 1 602 ] the  concen- 
tration of  oxygen- 16  in  the  permeating  gas,  [Si1802] 
the  concentration  of  oxygen-18  incorporated  in  the  silica 
scale,  and  a1  the  area  of  the  interstitial  sites. 

3.2.2  Case  I:  The  Rate  of  Permeation  Much  Greater  than 
the  Rate  of  Exchange 

It  the  rate  of  permeation  Is  much  greater  than  the 
rate  of  exchange  then  the  expected  oxygen-18  concentration 
versus  depth  profile  will  be  that  shown  in  Figure  15.  Here 


the  Incoming  oxygen-16  from  the  second  oxidation  has  a 
high  probability  of  being  transported  through  the  exist- 
ing silica  scale  and  reacting  at  the  silica-silicon 
Interface  without  exchanging  with  any  of  the  network 
oxygens.  It  Is  therefore  expected  that  the  oxygen-18 
will  remain  peaked  at  the  silica-oxygen  Interface. 

3.2.3  Case  II:  Rate  of  Permeation  Less  than  the  Rate 
of  Exchange 

In  this  case  the  thermal  oxidation  reaction  Is 
again  controlled  by  the  transport  of  molecular  oxygen  to 
the  silica-silicon  Interface,  but  the  relative  rates  of 
permeation  and  exchange  are  different  than  In  the  pre- 
vious case  causing  the  oxygen-18  marker  to  be  located  at 
the  silica-silicon  Interface.  This  Is  shown  In  Figure 
16.  The  development  of  this  profile  can  be  visualized  as 
follows:  as  a molecule  of  oxygen-16  dissolves  In  the 
silica  scale  the  probability  of  one  of  the  atoms  of 
oxygen-16  exchanging  with  an  oxygen-18  atom  in  the  net- 
work Is  greater  than  the  probability  of  It  moving  to  the 
next  site.  The  probability  of  the  other  oxygen-16 
exchanging  with  a network  oxygen-18  Is  also  high.  It  Is 
therefore  expected  that  the  oxygen-16  from  the  second 
part  of  the  oxidation  will  be  Incorporated  into  the 
silica  scale  at  the  silica-oxygen  Interface  and  the  oxygen 
reaching  the  silica-silicon  Interface  will  be  oxygen-18. 
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3.2.4  Case  III:  Relative  Rates  of  Permeation  and  Exchange 
Intermediate  to  Cases  I and  II 

It  can  be  expected  that  If  the  relative  rates  of 
exchange  and  permeation  are  Intermediate  to  Case  1 
(permeation  rate  much  greater  than  exchange  rate)  and 
Case  II  (permeation  rate  less  than  exchange  rate),  then 
the  predicted  oxygen-18  concentration  versus  depth  pro- 
file will  be  Intermediate  to  those  shown  in  Figures  15 
and  16.  Here  the  location  of  the  oxygen-18  "marker" 
will  be  peaked  between  the  silica-oxygen  Interface  and 
the  silica-silicon  Interface. 

The  development  of  this  type  of  profile  is  shown 
In  Figures  17a  - 17e  for  the  case  where  the  permeation 
rate  is  greater  than  the  exchange  rate  and  can  be  under- 
stood as  follows: 

Figure  17a  shows  the  expected  oxygen-18  concen- 
tration versus  depth  profile  at  the  completion  of  the 
Initial  oxidation.  Here  the  silica  scale  thickness  is 
y^  and  is  composed  entirely  of  oxygen-18. 

In  order  to  see  the  development  of  Figure  17b  - 
17e  it  Is  necessary  to  propose  some  type  of  model.  The 
simplest  model  Is  to  assume  that  the  silica  scale  can  be 
divided  into  discrete  layers  each  containing  the  same 
number  of  network  oxygens.  It  Is  then  assumed  that  the 
oxidation  reaction  occurs  In  discrete  steps  in  the  following 


Flour®  17t  Case  III:  The  relative  rates  of  permeatfor 
and  exchange  are  Intermediate  to  Cases  I and  II  at  various  oxid< 
thicknesses  (a)  y-y,;  (b)  y«1.25y,;  (c)  y*1.5y,;  (d)  y=1.75  y<; 
(e)  y*2. 0y4.  ' 1 1 1 


manner.  The  number  of  oxygens  required  to  form  the 
next  discrete  layer  of  oxide  Is  dissolved  In  the  first 
layer  of  silica.  This  assumes  an  oxygen  solubility  In 
silica  that  Is  much  greater  than  actually  exists  (the 
same  solution  Is  arrived  at  If  It  Is  assumed  a smaller 
amount  of  oxygen  Is  dissolved  In  the  first  layer  and 
the  process  repeated  until  enough  oxygen  has  passed 
through  the  layer  to  form  one  new  discrete  layer  of 
scale).  In  the  first  layer  a given  fraction  (depend- 
ing on  the  relative  rates  of  permeation  and  exchange) 
of  the  dissolved  oxygen  will  exchange  with  network 
oxygens.  This  exchange  will  occur  since  the  silica  is 
not  a perfect  network  of  silicon  and  oxygen,  l.e.  vacan- 
cies exist  in  the  silica  network.  These  vacancies  allow 
some  of  the  dissolved  oxygen  to  become  part  of  the  silica 
network,  and  to  maintain  thermodynamic  equilibrium  some 
of  the  network  oxygens  will  enter  the  dissolved  oxygen 
by  the  breaking  of  silicon  to  oxygen  bonds.  It  has  also 
been  assumed  here  that  there  is  no  Isotopic  preference 
to  the  exchange  reaction.  The  dissolved  oxygen  then 
moves  to  the  next  layer  where  the  process  repeats  itself. 
Once  the  dissolved  oxygen  has  traversed  the  entire  silica 
scale  it  reacts  completely  at  the  silica-silicon  inter- 
face to  form  a new  discrete  layer  of  scale.  The  amount 
of  oxygen  required  to  produce  the  next  layer  of  silica  is 
now  dissolved  in  the  first  layer  of  silica  and  the  above 


process  repeated.  This  continues  until  the  silica  scale 


thickness  has  Increased  to  2y^ . 

For  the  purpose  of  these  figures  It  will  be 
assumed  that  ten  percent  of  the  oxygen  exchanges  in  each 
layer.  Then  If  there  are  one  hundred  oxygens  In  each 
discrete  layer  we  find  that  one  tenth  of  the  oxygen  that 
dissolves  In  the  first  layer  to  form  the  new  oxide  will 
exchange.  Therefore  the  oxygen  moving  to  the  next  layer 
will  be  composed  of  ninety  atoms  of  oxygen-16  and  ten 
atoms  of  oxygen-18.  Ten  percent  of  these  oxygens  will 
now  exchange  with  the  second  layer  or  nine  oxygen-16's 
and  one  oxygen-18.  The  gas  moving  to  the  third  layer  is 
then  eighty-one  percent  oxygen-16  and  nineteen  percent 
oxygen-18.  This  process  continues  until  one-tenth  of 
the  dissolved  oxygen  has  exchanged  with  the  oxygen  in  the 
last  layer  of  silica.  This  oxygen  then  reacts  with  the 
silicon  to  form  new  oxide  with  a composition  equal  to  the 
composition  of  the  dissolved  gas  after  the  final  exchange. 
Once  this  layer  of  oxide  Is  formed  the  oxygen  required  to 
form  the  second  layer  of  new  oxide  will  be  dissolved  in 
the  first  layer.  Here  one-tenth  or  ten  of  the  oxygens  will 
exchange  with  the  network  oxygens,  but  the  first  layer 
presently  has  a composition  of  ten  percent  oxygen-16  and 
ninety  percent  oxygen-18  making  only  ninety  percent  of 
these  exchanges  effective.  After  the  exchanges  occur  the 
scale  will  be  eighty-one  percent  oxygen-18  and  the  gas 
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moving  to  the  second  layer  will  be  nine  percent  oxygen-18. 
Ten  percent  of  the  oxygen-16  in  the  gas  will  now  exchange 
with  the  oxygen  in  the  second  layer,  but  only  ninety-one 
percent  of  these  exchanges  will  effect  the  scale  composi- 
tion. Ten  percent  of  oxygen-18  in  the  gas  will  also 
exchange  with  the  silica  scale  and  nine  percent  of  these 
exchanges  will  be  effective  in  changing  the  scale  com- 
position. The  new  scale  will  now  be  composed  of  eighty- 
four  percent  oxygen-18  and  sixteen  percent  oxygen-16.  The 
dissolved  gas  transferred  to  the  third  layer  will  be  six- 
teen percent  oxygen-18.  This  process  continues  until  the 
permeating  gas  has  exchanged  with  the  last  silica  layer. 

It  then  reacts  with  the  silicon  forming  a new  layer  of 
silica  and  the  gas  required  for  the  next  silica  layer  is 
dissolved  and  the  process  repeated  until  the  scale  thick- 
ness has  Increased  to  2y^ . 

Figures  17b  - e show  the  calculated  oxygen-18  con- 
centration versus  depth  profiles  assuming  that  ten  percent 
of  the  oxygen  exchanges  In  each  layer.  Figure  17b  shows 
the  oxygen-18  concentration  versus  depth  profile  when  the 
scale  thickness  has  increased  to  1.25  y..  It  can  be  seen 
the  oxygen-18  concentration  exhibits  a maximum  within  the 
scale  Instead  of  at  either  interface.  Figure  17c  shows 
the  concentration  versus  depth  profile  when  the  scale  is 
1.5  y^  thick.  Figures  17d  and  17e  show  the  continued 


development  of  the  oxygen-18  concentration  versus  depth 
profile  when  the  scale  thickness  has  Increased  to  1.75  y,. 
and  2.0  y^  respectively.  It  can  be  seen  from  Figure  17e 
that  the  final  profile  exhibits  a maximum  In  oxygen-18 
concentration  profile  at  approximately  1.3  y^. 

3.3  Oxidation  Controlled  by  the  Transport  of  an  Ionic  Species 

3.3.1  Case  IV:  Ionic  Transport  of  Oxygen 

If  the  thermal  oxidation  of  silicon  is  controlled 
by  the  transport  of  an  Ionized  species  of  oxygen  as  pre- 
dicted by  Jorgensen  (13)  the  concentration  profile  will  be 
given  by  the  Wagner  Theory  (see  Section  2.1.4)  or  the  Morin 
modification  to  the  theory  of  point  defect  diffusion  (52). 
Here  the  oxygen-18  will  be  found  peaked  at  the  silicon- 
silica  interface  as  Is  shown  in  Figure  18. 

3.3.2  Case  V:  Ionic  Transport  of  Silicon 

It  Is  also  possible  that  the  thermal  oxidation 
reaction  Is  controlled  by  the  diffusion  of  silicon  ions 
with  the  reaction  occurring  at  the  si  1 ica* oxygen  inter- 
face. If  this  is  the  situation  the  oxygen-18  from  the 
first  oxidation  is  expected  to  be  located  at  the  silicon- 
silica  Interface  as  a step  function  to  first  approxi- 
mation. The  expected  oxygen-18  distribution  is  shown  in 
Figure  19. 
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EXPERIMENTAL  PROCEDURE 

4 . 1 Samples 

The  silicon  coupons*used  in  these  thermal 
oxidation  studies  were  oriented  in  the  (110)  direction 
and  were  n-type  crystals  with  a phosphorus  dopant  level 

1 5 

of  10  atoms  per  cubic  centimeter.  The  coupons  were 
one  half  millimeter  thick  and  cut  to  an  eight  millimeter 
square. 

The  coupons  were  mirror  finished  on  one  side  and 
were  cleaned  by  the  following  procedure:  (1)  deionized 
water  rinse;  (2)  ultrasonic  rinse  in  NH40H:  H202:  H20 

(1:1:5);  (3)  repeat  1;  (4)  ultrasonic  rinse  in  HC1:H202: 
H20  (1:1:5);  (5)  repeat  1;  (6)  48%  HF  rinse;  (7)  repeat 
1;  (8)  dry  in  dried  air. 

4.2  Thermal  Oxidations 

4.2.1  Oxidation  Apparatus  and  General  Procedures 

Figure  20  shows  a schematic  diagram  of  the 
oxidation  apparatus.  The  system  is  designed  in  such 


* This  single  crystal  silicon  was  obtained  from  Ventron 
Corporation  of  Bradford,  Pennsylvania. 
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Figure  20:  Schematic  diagram  of  the  oxidation  apparatus. 


a way  as  to  allow  the  static  oxidation  of  up  to  five 
samples  which  rested  on  alumina  shelves  In  the  sample 
holder.  A Pt-6%Rh,  Pt-30%  Rh  thermocouple  Is  located 
next  to  each  sample  position  to  enable  continuous 
temperature  monitoring. 

The  oxygen  used  In  the  oxidation  experiments  was 
stored  In  Pyrex  bulbs  which  have  a zeolite-filled  finger 
to  aid  In  removing  the  oxygen  from  the  system  and  to  In- 
sure that  the  oxygen  Is  dry.  The  pure  oxygen  gas  was  en- 
riched to  between  10  and  20  atom  percent  oxygen-18;  and 
the  oxygen-18  depleted  gas  was  research-grade  oxygen. 

After  placing  the  silicon  coupons  In  the  oxidation 

chamber,  the  chamber  was  evacuated  using  a liquid-nitrogen- 

trapped  mercury  diffusion  pump  to  a pressure  of  no  more 

than  10  microns  of  mercury.  The  system  was  then  leak 

checked  using  a helium  leak  detector  to  Insure  that  the 

system  was  leak  free.  The  oxidation  chamber  was  then 

heated  to  temperature  while  continued  pumping  maintaining 

★ 

a vacuum  of  no  more  than  10  microns  of  mercury.  Once 
the  system  reached  the  desired  temperature  it  was  again 
leak  checked  to  be  sure  no  leaks  had  developed  during 
heating.  If  no  leaks  were  present,  the  system  was  pres- 
surized to  approximately  780  millimeters  of  mercury  with 
oxygen.  During  the  oxidation  run,  the  sample  temperature 


was  contlnously  monltered  using  a strip  chart  recorder. 
When  the  desired  oxidation  time  was  reached,  the  furnace 
was  shut  off  and  the  oxygen  removed  via  the  zeolite 
fingers.  Once  the  oxygen  had  been  removed  the  sample 
holder  was  raised  and  the  furnace  lowered  to  allow  the 
system  to  finish  cooling.  When  the  oxidation  chamber 
reached  a temperature  of  approximately  800*K  dry  air  was 
Introduced  to  the  chamber.  Thi:  entire  shut-down  pro- 
cedure required  thirty  minutes  from  the  time  the  furnace 
was  shut  off  until  the  samples  could  be  transferred  to 
a desslcator. 

Throughout  the  course  of  the  oxidation  run,  gas 
samples  were  periodically  withdrawn  in  order  to  monitor 
the  Isotopic  composition.  A minimum  of  four  gas  samples 
were  drawn  during  each  run. 

4.2.2  Kinetic  Oxidation  Experiments 

The  kinetic  experiments  involved  oxidizing  sam- 
ples In  an  oxygen-18  enriched  oxygen  atmosphere  for 
various  times  at  several  temperatures.  The  temperature 
range  selected  for  the  kinetic  study  was  1089±5°K  to 
1473±5°K.  The  oxidation  times  ranged  from  ten  to  thirty 
hours . 

4.2.3  Tracer  Experiments 

The  tracer  experiments  were  preformed  in  the  same 


temperature  range  as  the  kinetic  experiments.  The 
silicon  coupons  were  oxidized  In  an  oxygen-18  enriched 
oxygen  atmosphere  to  give  a silica  scale  thickness  of 
approximately  4500  A.  The  oxidation  was  the  Interrupted 
and  the  oxygen-18  enriched  oxygen  replaced  with  oxygen- 
18  depleted  oxygen  and  the  oxidation  continued  to  approx- 
imately 9000  A. 

These  experiments  were  repeated  with  the  isotopic 
order  reversed  as  an  Internal  check  on  the  experiment. 

4.3  Charged  Particle  Activation  Analysis 

4.3.1  Sample  Preparation  and  Positioning 

After  the  coupons  were  oxidized  they  were  placed 

in  the  center  of  an  aluminum  blank  one  inch  in  diameter 

and  of  a thickness  such  that  the  combined  coupon-blank 

thickness  was  0.275  inch.  This  assemblage  was  then 

placed  in  an  evaporator  and  a thin  layer  of  aluminum 

* 

(approximately  100  A)  deposited  over  the  sample  surface. 
The  samples  were  then  sandwiched  between  two  aluminum 
blanks  with  the  top  blank  having  a 0.25  inch  diameter 
hole  in  the  center  beveled  to  0.5  inch  diameter  at  the 
top  of  the  blank. 


* The  aluminum  film  thickness  was  calculated  assuming  4„ 
distribution  from  the  tungsten  holder. 


This  sample  preparation  has  proven  to  be  adequate 
in  conducting  both  heat  and  charge  from  the  sample  sur- 
face during  proton  bombardment. 

Once  the  samples  have  been  secured  in  their  alu- 
minum holders,  they  were  mounted  in  a movable  frame, 
capable  of  holding  five  samples,  in  the  scattering  cham- 
ber. The  frame  Is  positioned  in  such  a way  that  the 
normal  to  the  sample  surface  makes  an  angle  of  7»s  degrees 
with  both  the  incident  proton  beam  and  the  alpha  particles 
(the  tilting  of  the  sample  was  not  necessary  in  these 
studies,  but  is  of  importance  when  working  with  single 
crystals  to  prevent  channelling  of  the  protons  or  alphas). 
The  scattering  chamber  was  then  sealed  and  a vacuum  drawn. 
The  chamber  was  maintained  at  a pressure  of  no  more  than 
70  millimicrons  of  mercury  throughout  the  analysis. 

4.3.2  Alpha  Particle  Measurement  and  Recording 

The  energetic  alpha  particles  observed  at  a 
laboratory  scattering  angle  of  165°  pass  through  a mag- 
netic spectrometer,  designed  by  Lindstrom  and  Heuer  (46), 
which  removes  the  backscattered  protons  from  the  alpha 
particles  by  a magnetic  field.  The  magnetic  spectro- 
meter eliminates  the  need  for  mylar  to  separate  the 
alphas  from  the  backscattered  protons,  which  is  an  advan- 
tage since  one  does  not  have  to  worry  about  the  uniformity 
or  stopping  power  of  the  mylar  (the  magnetic  field  only 


changes  the  particle's  direction  and  does  not  alter 
their  energy).  Figure  21  shows  a schematic  of  the 
scattering  chamber  and  magnetic  spectrometer  used 
in  these  experiments. 

A gold-coated  surface  barrier  detector  was 
used  to  record  the  energetic  particles.  The  detector 
signal  was  processed  (amplified  and  shaped)  and  sub- 
sequently recorded  in  a multichannel  pulse-height 
analyzer  using  a standard  electronic  set-up  (Figure  22). 

At  the  conclusion  of  the  run  each  recorded 
spectrum  was  transferred  to  punched  paper  tape  for 
later  computer  analysis. 

4.3.3  The  Activation  Experiment 

The  multichannel  analyzer  was  calibrated  at 

1 8 

the  beginning  and  end  of  each  day  using  a Ta2  05  target* 

**  18 

and  a ZnF2  target.  The  Ta2  05  calibration  run  was  made 
at  an  incident  proton  energy  of  750  KeV  and  the  ZnF2  run 
at  850  KeV  and  890  KeV  giving  three  calibration  points. 

Once  the  MCA  had  been  calibrated  the  samples 
were  bombarded  with  protons,  and  the  alpha  spectrum 

* The  Ta21B0s  targets  were  prepared  by  Dr.  J.P.S.  Pringle,  Chalk 
River  Nuclear  Laboratories,  Atomic  Energy  of  Canada,  Ltd.,  Chalk 
River,  Ontario,  Canada,  by  anodically  oxidizing  tantalum  blanks. 

**The  ZnF2  targets  were  prepared  by  evaporating  anhydrous  ZnF2 
onto  tantalum  blanks. 
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Schematic  of  the  scattering  chamber  and  magnetic  spectrometer. 
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Figure  22:  Schematic  of  elctronic  circuitry 
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recorded  for  a predetermined  proton  flux.  A Ta2  05 
sample  was  run  between  each  sample  run  in  order  to 
determine  the  stability  of  the  detector  and  current 
integrator.  Each  sample  was  run  at  an  incident  pro- 
ton energy  of  750  KeV  for  an  uncorrected  integrated 
charge  of  1000  microcoulombs,  and  at  an  incident 

energy  of  800  KeV  for  an  integrated  charge  of  500 

18 

microcoulombs.  The  interspersed  Ta2  05  runs  were  100 
microcoulombs  duration  and  at  the  same  incident  energy 
as  the  proceding  run.  A beam  current  of  0.5  micro- 
amphere  was  used  for  all  runs. 

The  linearity  of  the  MCA  and  electronics  was 
periodically  checked  using  a pulse  generator.  The 
connection  of  the  pulse  generator  to  the  system  is 
shown  In  Figure  22. 

4.3.4  Current  Integrator  Calibration 

The  current  integrator  system  was  also  per- 
iodically checked  by  placing  a known  current  across 
the  target  and  measuring  the  time  required  to  collect 
a given  amount  of  charge.  The  system  used  to  calibrate 
the  current  integrator  is  shown  schematically  in  Figure 
23. 
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SECTION  V 


I 


DATA  ANALYSIS 

5 . 1 Current  Integrator  Cal i brat  ion 

Periodic  calibration  of  the  current  integrator 
(see  Section  4.3.4)  demo'strated  excellent  day-to-day 
reproducibility.  The  true  integrated  charge,  Qtrue 
was  found  to  be  related  to  that  observed,  QQbs  by 

Qtrue  = 1*01"(Qobs  " °-003  (42) 

where  t is  the  time  in  seconds  required  to  collect  Qobs. 
The  multiplicative  factor  is  the  calibration  and  the 
negative  correction  term  the  offset  of  the  charge  pump 
circuit  used. 

5.2  MCA  Cal i brati on 

1 8 

5.2.1  Ta2  O5  Calibration 

1 8 

Each  of  the  alpha-spectra  from  the  Ta2  05  sample 

was  initially  fitted  to  a Gaussian  function.  This  was 

possible  since  the  unconvoluted  spectra  obtained  from  a 
> 1 8 

530  A thick  Ta2  05  sample  would  be  approxi ma tely  15  KeV 

thick  in  alpha  energy  (this  can  be  seen  from  the  range- 

energy  curves  for  protons  and  alphas  on  Ta205  and  the 

1 8 , , 1 s 

kinematics  of  the  0(p,a)  N reaction  given  in  Appendix 
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1).  The  channel  corresponding  to  the  center  of  the 
Gaussian  peak  was  therefore  the  channel  which  corre- 
sponded to  the  energy  of  an  alpha  particle  created 

it) 

at  the  center  of  the  Ta2  05  sample.  The  energy  of 

this  alpha  particle  Is  easily  calculated  from  the 

18 

range-energy  curves  for  protons  and  alphas  on  Td2  05, 

18  is 

the  kinematics  of  the  0(p,a)  N reaction,  and  the 
half  thickness  of  the  oxide  sample. 

5.2.2  ZnF2  Calibration 

Each  ZnF2  spectrum  was  fitted  In  the  same 
18 

way  as  the  Ta2  05  spectra.  Protons  of  Incident 
energy  approximately  850  KeV  undergo  the  nuclear 

19  . .16  * 

reaction  F(p,oj)  0 . The  energetic  alpha  particle 
corresponding  to  the  ground  state  of  oxygen-16  has 
approximately  6.66MeV  of  energy  and  that  correspond- 
ing to  the  first  excited  state  of  oxygen-16  has 
approximately  1.99  MeV  of  energy. 

The  range-energy  curves  for  protons  and  alphas 

19  16 

on  ZnF2  and  the  kinematics  of  the  F(p,a)  0 and 

19  , 16  * 

F(p,cn)  0 reactions  are  given  In  Appendix  1. 

5.2.3  MCA  Calibration  Channel,  MCA  Dispersion,  and 
Spreading  Function 

1 8 

The  alpha  energies  from  the  Ta2  05  and  ZnF2 
calibrations  were  graphed  as  a function  of  channel 
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number.  As  can  be  seen  from  Figure  24  the  alpha  energy 
Is  a linear  function  of  channel  number.  The  MCA  cali- 
bration then  consisted  of  one  point  from  the  graph  and 
the  slope  and  dispersion  In  KeV  per  channel. 

The  calibration  channel  and  dispersion  deter- 

1 8 

mined  above  was  then  used  to  fit  several  of  the  Ta2  05 

alpha  spectra  to  obtain  the  spreading  function.*  The 

correct  spreading  function  was  found  by  using  a non- 

2 

linear  regression  scheme  based  on  the  x test  (53), 

2 

the  minimum  reduced  chi  square  , xr»  giving  the  best 
value  of  the  spreading  function. 

5.2.4  Pulser  Calibration 

Each  of  the  pulser  peaks  obtained  from  the 

pulser  calibration  was  fitted  to  a Gaussian  function  in 

16 

the  same  manner  as  the  Ta2  05  and  ZnF2  calibration  peaks 
-(The  pulser  peak  is  spread  only  by  the  electronics  of 
the  system).  The  pulser  setting  was  then  plotted  as  a 
function  of  channel  number  to  check  the  linearity  of  MCA. 
A typical  series  of  pulser  runs  Is  shown  In  Figure  25. 

5 . 3 Fitting  the  Alpha  Spectra 
5.3.1  General 

The  recorded  alpha  spectrum  for  each  of  the  sam- 


* The  computer  program  was  written  by  Dr.  Walter  W.  Llndstrom. 
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Channel  430.80  is  3359.66  KeV  with 
6.33  KeV  per  channel 


19  , 16  * 

F(p,ai)  0 


MCA  Channel 

Figure  24:  Graph  used  to  calibrate  MCA 


pies  was  fitted  knowing  the  KCA  calibration  (determined 

above)*  the  spreading  function,  and  the  range-energy 

curves  for  protons  and  alphas  on  amorphous  silica  (these 

curves  are  given  In  Appendix  1)  using  a trial  and  error 

convolution  procedure.  In  the  kinetic  experiments  only 

the  thickness  of  the  silica  layer  was  varied  and  In  the 

mechanistic  experiments  only  the  concentration  versus 

depth  profile  was  varied  to  fit  the  experimental  alpha 

spectra.*  The  experimental  spectra  were  fit  using  a 

2 

non-linear  regression  scheme  based  on  the  x test  (53), 

2 

with  the  best  fit  being  given  by  the  minimum  x- 

■ • 

5.3.2  Kinetic  Oxidation  Fits 

It  was  possible  to  fit  the  alpha  spectra  from 
the  kinetic  oxidation  experiments  with  a step  function 
oxygen-18  profile,  since  the  oxygen-18  concentration  In 
the  ambient  gas  remained  constant  throughout  the  runs. 

Each  of  the  recorded  spectrum  was  therfore  fit  assuming 
a constant  oxygen-18  concentration  within  the  silica 
scale.  This  concentration  was  found  nearly  equal  to  the 
gas  concentration.  Implying  a very  rapid  surface  exchange. 

5.3.3  Mechanistic  Oxidation  fits 

The  mechanistic  oxidation  experimental  data  was 

fitted  using  a trial  and  error  convolution  procedure. 

* The  computer  program  was  written  by  Dr.  Walter  W.  Llndstrom. 
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In  each  case  It  was  Initially  assumed  that  the  sample 

had  a uniform  oxygen-18  distribution  and  the  predicted 

spectrum  compared  with  the  experimental  data;  the  oxygen- 

18  distribution  was  then  modified  and  the  procedure 

2 

repeated  until  a reasonable  xr  was  obtained.  The  pro- 
cedure used  to  fit  the  data  from  one  of  the  samples  Is 
shown  In  Appendix  2. 

Typical  experimental  spectra  and  the  correspond- 
ing calculated  spectra  are  shown  In  Appendix  3 for  some  of 
the  mechanistic  experiments. 


SECTION  VI 

RESULTS  AND  DISCUSSION 
6.1  Kinetic  Experiments 

The  square  of  the  measured  silica  scale  thickness, 
as  measured  by  the  single  spectrum  proton  activation  tech- 
nique, was  plotted  as  a function  of  time  for  each  of  the 
temperatures  Investigated  and  the  parabolic  rate  constant 
determined  from  the  slope  of  the  line  (Figures  26-29). 

Note  that  the  point  (0,0)  was  used  In  determining  these 
rate  constants  since  at  time  t * 0 there  Is  no  oxygen-18 
silica  scale.  The  logarithm  of  these  rate  constants  was 
then  plotted  as  a function  of  reciprocal  temperature  to 
determine  the  general  expression  for  the  parabolic  rate 
constant  In  the  temperature  range  1089±5°K  to  1473±5PK. 

This  graph  Is  shown  In  Figure  30  and  a least  squares  fit 
gives  the  following  general  expression  for  the  parabolic 
rate  constant: 

kp  * [(7.5±3.0)  x 10"10  exp  - (^8ffi°-9eV)]  (42) 

The  experimentally  determined  activation  energy 
found  In  these  kinetic  studies,  1.08±0.09eV,  Is  0.13eV 
lower  than  any  value  reported  In  the  literature  for  the 
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t * 10  (sec) 


Figure  26:  The  measured  silica  scale 
Thickness  squared  as  a function  of  time  at 
1091  ± 5°K. 


10"  (cm  ) 
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of  silicon  In  "dry"  oxygen.  (The 
maximum  value  of  the  activation  energy  Is  one  half  standard 
deviation  lower  than  the  lowest  value  reported  In  the  liter- 
ature). The  silica  scale  thicknesses  as  measured  In  these 
experiments  were  found  to  be  consistently  thicker,  by  ten 
to  twenty  percent,  than  other  reported  measurements  for  the 
thermal  oxidation  of  silicon.  However,  It  was  found  that 
the  scale  thickness  measured  by  this  technique  agreed  with 
the  estimated  thickness  corresponding  to  the  Interference 
color  of  the  scale.  The  correlation  of  the  measured  thick- 
ness and  estimated  thickness,  using  the  color  chart  given  by 
Pliskln  and  Conrad  (54)  based  on  the  VAMFO  * technique  for 
measuring  thin  films.  Is  given  In  Table  4. 

The  reduced  chi  square  values  obtained  for  the 
kinetic  sample  fits  used  In  the  above  calculations  ranged 
from  1.3  to  3.0  and  were  based  on  only  the  MCA  channels 
which  contained  non-zero  data  points.  (The  Inclusion  of  MCA 
channels  which  are  expected  to  have  zero  data  points  leads 
to  better  reduced  chi  square  values;  however,  these  are 
meaningless  In  determining  the  goodness  of  the  fit.  For 
example  a spectrum  containing  forty  non-zero  data  points 
may  have  a reduced  chi  square  of  1.86,  but  If  one  Includes 
70  data  points  of  which  30  are  zero's  the  reduced  chi  square 


* VAMFO  stands  for  Variable  Angle  Monochromatic  Fringe  Observation. 


TABLE  4.  —Correlation  of  the  measured  thickness  and  the  estimated  thick 
ness,  using  the  color  chart  given  by  Pllskln  and  Conrad  (54). 


S- 11-2  7.2x10  1095  0.24  Yellow-Green  0.17-0.22 


becomes  1.01).  Figure  31  shows  the  computer  predicted  fit 
and  data  points  for  sample  S-10-1  (0.66  micron  thick) 
which  was  oxidized  for  3.6  * 10  seconds  (10  hours)  at 
1337±5#K  and  analyzed  using  protons  of  Incident  energy 
750  KeV.  The  calculated  reduced  chi  square  value  for  this 
fit  is  1.86.  It  should  be  noted  that  the  region  of  the 
spectra  most  difficult  to  fit  is  the  region  correspondi ng 
to  the  center  of  the  scale.  The  program  consistently 
predicted  too  low  an  alpha  yield.  The  fact  that  the  pre- 
dicted alpha  yield  at  the  center  of  the  scale  is  too  low 
is  an  Indication  the  stopping  powers  for  protons  and  alphas 
assumed  for  silica  may  be  in  error.  The  effect  of  varying 
the  proton  and  alpha  stopping  powers  will  be  treated  below. 

The  lower  activation  energy  found  in  these  experi- 
ments as  well  as  the  thicker  thermal  oxide  scales  can  pos- 
sibly be  attributed  to  one  or  more  of  the  following:  (1) 
the  water  content  of  the  oxygen,  (2)  the  variation  of  the 
thermal  oxide  density  with  temperature  of  formation,  and 
(3)  an  error  in  the  stopping  power  for  protons  and  alphas 
assumed  for  silica.  Each  of  these  will  be  discussed  in 
turn. 

It  has  been  reported  by  Irene  (17)  that  the  para- 
bolic activation  energy  can  be  varied  by  as  much  as  0.2  eV 
by  changing  the  water  vapor  content  of  the  ambient  from 
less  than  one  part  per  million  to  twenty-five  parts  per 
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Figure  31:  Fit  of  sample  S-10-1,  0.66y  thick,  at  an  Incident  proton  energy  of  750  KeV 


million.  The  oxygen  Introduced  to  the  oxidation  system 
has  a dew  point  of  208°K  which  corresponds  to  a water  vapor 
content  of  approximately  thirteen  parts  per  million.  The 
water  vapor  content  of  the  ambient  oxygen  can  also  be  expect- 
ed to  Increase  during  the  course  of  the  run  due  to  the  oxida- 
tion of  any  organic  molecules  In  the  system.  Trace  amounts 
of  carbon  dioxide,  also  a product  of  organic  oxidation, 
have  been  found  In  the  gas  samples,  Indicating  this  to  be  « 
possible  source  of  additional  water  vapor  In  the  ambient. 
Although  the  dew  point  of  the  oxygen  was  not  determined 
during  the  experiment,  an  estimate  of  the  maximum  water 
vapor  content,  based  on  the  equilibrium  values  of  the 
zeolite-water  system.  Is  30  parts  per  million.  Never- 
theless, this  amount  of  water  vapor  could  effect  the 
measured  activation  energy  by  as  much  as  three  percent. 

It  was  assumed  In  these  studies  that  the  density  of 
the  thern.al  silica  scale  was  2.2C  grams  per  cubic  centimeter. 
However,  Deal  (15)  has  found  that  the  density  of  the  thermal 
silica  layer  formed  at  1273°K  Is  2.27  grams  per  cubic  centi- 
meter and  at  1473°K  Is  2.15  grams  per  cubic  centimeter. 

Since  the  stopping  power  of  a particle  In  matter  Is  expressed 
In  units  of  energy-area  per  unit  mass,  a variation  in  thermal 
oxide  density  would  therefore  effect  the  range-energy  curve 
for  the  particle.  This  variation  In  density  could  introduce 
a systematic  error  in  the  scale  thickness  of  as  much  as  four 
percent  at  each  temperature. 
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The  elemental  stopping  powers  used  to  calculate  the 
stopping  power  of  silica  Is  a third  possible  source  of 
systematic  error.  The  elemental  stopping  powers  as  given 
by  Williamson,  Boujot,  and  Picard  (55)  can  be  considered 
accurate  to  no  more  than  ten  percent.  An  Increase  In  pro- 
ton stopping  power  would  tend  to  shift  the  energy  of  alphas 
created  below  the  surface  to  lower  energies  while  an  Increase 
In  alpha  stopping  power  would  tend  to  shift  the  energy  of 
alphas  created  below  the  surface  to  lower  energies.  It 
should  also  be  noted  that  changing  the  proton  stopping  powers 
will  not  effect  the  general  shape  of  the  spectra,  but  will 
change  the  total  alpha  yield.  This  Is  because  the  only  term 
effected  In  the  equation  for  calculating  the  alpha  yield 
(Equation  37)  Is  the  differential  cross  section  In  a given 
depth  Interval.  A change  In  alpha  stopping  power  will  not 
effect  the  overall  alpha  yield  (this  assumes  that  all  alpha 
particles  created  have  enough  energy  to  reach  the  detector 
and  be  recorded  in  the  MCA),  but  will  change  the  spectral 
shape  due  to  the  rate  at  which  the  alpha  particles  lose 
energy  as  they  traverse  the  sample.  Figures  32-37  show  the 
effect  of  varying  the  proton  and  alpha  stopping  powers,  by 
plus  or  minus  ten  percent  either  singularly  or  simultane- 
ously, on  the  calculated  alpha  spectra  for  a silica  scale 
of  uniform  oxygen-18  concentration  0.66  micron  thick. 

It  can  be  seen  from  Figures  32-37  that  the  predicted 
alpha  yield  at  the  center  of  the  scale  can  be  increased  by 
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gure  32:  The  effect  of  changing  the  proton  stopping 


Figure  35:  The  effect  of  Increasing  the  proton  stopping  power  by  10*  while 
decreasing  the  alpha  stopping  power  by  10*. 


Channel  Number 

Figure  36:  The  effect  of  decreasing  both  the  proton  and  alpha  stopping  powers 


Figure  37:  The  effect  of  decreasing  the  proton  stopping 
Increasing  the  alpha  stopping  power  by  10*. 


1 


one  of  the  following:  (1)  Increasing  the  proton  stopping 
power,  (2)  decrea**n9  the  a1Pha  stopping  power,  (3)  Increas- 
ing the  proton  st"f'Pln9  power  while  decreasing  the  alpha 
stopping  power,  o»  (*)  decreasing  both  the  proton  and  alpha 
stopping  powers.  ln  or<ler  to  determine  how  the  stopping 
powers  should  be  v-rled  to  Improve  the  spectral  fits  It  Is 
necessary  to  have  6 series  of  silica  scaies  of  known  thick- 
ness and  of  unlfo*1"  oxygen-18  distribution.  These  s*.vles 
should  then  be  ru"  at  several  Incident  proton  energies 
where  there  are  l^ge  variations  In  cross  section  with 
Incident  proton  (such  as  the  637  KeV  resonance). 

This  will  give  a ceries  of  alpha  spectra,  each  having  a dif- 
ferent alpha  yiel'1  and  different  spectral  shape.  Since  the 
thickness,  oxygen  concentration,  and  cross  section  are 
known  It  will  be  r*»s s 1 bl e to  fit  these  spectra  by  varying 
only  the  proton  a*"*  *1Pha  stopping  powers.  Computer  calcula- 
tion shows  that  a variation  of  up  to  ten  percent  In  proton 
and/or  alpha  stopi,in9  Powers  can  cause  an  error  of  approxi- 
mately ten  percen*  I*1  the  measured  scale  thickness. 

6.2  "Marker"  Expe!illie-nts 

In  the  "marker"  experiments  It  was  found  that  the 
position  of  the  ,.-voen-18  "marker"  varied  with  oxidation 
temperature.  In  « he  temperature  range  1171±5°K  to  1473±5°K 


the  position  of  Isotope  used  In  the  Initial  oxidation 


moved  from  the  silica-oxygen  Interface  to  the  silicon- 
silica  Interface  with  Increasing  oxidation  temperature. 

The  oxygen-18  concentration  versus  depth  profiles  were 
derived  using  a trial  and  error  deconvolution  procedure. 

Each  spectrum  was  Initially  fit  assuming  a silica  scale 
of  uniform  oxygen-18  concentration  (20%)  and  one  micron 
thick.  The  spectra  correspond  .ig  to  this  oxygen-18  dis- 
tribution was  then  compared  with  the  experimental  data. 
Figure  38  shows  a comparison  of  the  experimental  data  for 
sample  S-5-3  analyzed  at  750  KeV  and  predicted  alpha  spectra 
assuming  a uniform  oxygen-18  distribution. 

Knowing  the  Incident  proton  energy,  the  kinematics 
18  15 

of  the  0(p,a)  N reaction  (Appendix  1),  the  range-energy 
curves  for  protons  and  alphas  on  silica  (Appendix  1),  and 
the  MCA  calibration,  it  Is  a simple  procedure  to  determine 
what  depth  each  channel  of  the  MCA  corresponds  to.  If  the 
experimental  and  predicted  spectra  were  Ideal  (l.e.  no 
spreading  occurred)  then  the  number  of  counts  In  a given 
channel  of  the  experimental  data  divided  by  the  number  of 
counts  In  the  same  channel  of  the  assumed  spectra  times  the 
concentration  of  oxygen-18  assumed  to  be  at  the  depth  corre- 
sponding to  the  channel  would  be  the  concentration  at  that 
depth.  Since  the  spectra  are  spread,  this  value  can  only 
be  used  as  a first  approximation  to  the  oxygen-18  concen- 
tration versus  depth  profile.  With  this  approximation  it  was 
possible  to  obtain  an  Improved  concentration  versus  depth 


Figure  38:  Initial  fit  of  sample  S-5-3  at  an  Incident  proton  energy  of  750  KeV. 


profile.  Continued  Iteration  of  this  procedure  allowed  the 
predicted  spectra  to  assume  the  same  shape  as  the  experi- 
mental data.  Figure  39  shows  the  final  fit  for  sample  S-5-3 
at  750  KeV.  (Appendix  2 Illustrates  the  entire  procedure 
used  to  fit  sample  S-5-3  at  both  750  Kev  and  800  KeV).  Using 
this  procedure  it  was  possible  to  obtain  reduced  chi  square 
values  which  ranged  from  1.65  to  4.0.  When  the  reduced  chi 
square  values  reached  this  range  each  Iteration  yielded  only 
a small  Improvement  In  the  reduced  chi  square  value.  At 
this  point  farther  refinement  of  the  concentration  versus 
depth  profile  became  dependent  on  the  amount  of  time  avail- 
able for  calculating  each  new  assumed  profile.  To  obtain 
reduced  chi  square  values  less  than  2.0  In  a reasonable 
amount  of  time  would  require  a deconvolution  program. 

It  Is  also  possible  to  Infer  something  about  the 
location  of  the  oxygen-18  "marker"  by  looking  at  the  experi- 
mental alpha  spectra.  This  can  be  seen  by  comparison  of  the 
experimental  data  for  two  samples  oxidized  under  Identical 
conditions,  with  the  exception  of  the  Isotopic  order  of 
oxidation,  and  analyzed  at  the  same  Incident  proton  energy. 
Figure  40  shows  the  experimental  data  for  sample  S-4-1 
oxidized  In  oxygen-18  followed  by  oxygen-16  at  1473±5°K  and 
analyzed  at  an  Incident  energy  of  750.0  KeV.  Figure  41 
shows  the  experimental  data  for  sample  S-5-1  treated  Identically 
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Figure  41:  The  experimental  alpha  spectrum  for  sample  S-5-1  at  an  Incident 
energy  of  750  KeV. 


to  sample  S-4-1  with  the  exception  of  the  isotopic 
order.  By  looking  at  the  front  edge*  high  energy  side.  It 
can  readily  be  seen  that  the  data  points  In  Figure  41, 
sample  S-5-1,  rise  more  sharply  than  do  those  In  Figure  40, 
sample  S-4-1.  This  Indicates  that  the  oxygen-18  In  sample 
S-5-1  Is  located  closer  to  the  surface  than  Is  the  oxygen- 
18  In  sample  S-4-1.  Similarly  by  looking  at  the  trailing 
edge,  low  energy  side,  it  can  be  seen  that  the  oxygen-18 
In  sample  S-4-1  Is  concentrated  deeper  In  the  sample  than 
for  sampl e S-5-1 . 

There  are  two  possible  sources  of  error  In  determin- 
ing the  concentration  versus  depth  profiles  for  these  samples. 
They  are  (1)  the  variation  In  the  density  of  the  thermally 
formed  silica  with  temperature  and  (2)  an  error  In  the  stop- 
ping power  for  protons  and  alphas  In  silica.  A variation  in 
the  density  of  the  thermally  formed  silica  will  effect  only 
the  apparent  thicknesr  of  the  thermal  scale,  but  will  not 
effect  the  relative  location  of  the  "marker".  The  same  Is 
true  If  there  Is  an  error  In  the  stopping  powers  assumed 
for  protons  and  alphas  in  silica.  This  can  be  seen  in 
Figures  32-37.  These  two  errors  will  therfore  only  effect 
the  goodness  of  the  fit  and  not  the  relative  "marker"  loca- 
tion that  the  results  are  based  on. 

Figures  42-45  show  the  fractional  oxygen-18  con- 
centration versus  depth  profiles  for  the  samples  that  were 
oxidized  In  oxygen-18  followed  by  oxygen-16.  The  concentration 
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Figure  43:  Concentration  versus  depth 
profile  for  sample  S-4-2  oxidized  at  1412±5°K. 
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versus  depth  profiles  for  the  samples  where  the  Isotopic 
order  was  reversed  are  not  shown  since  they  are  simply 
Identical  profiles  with  the  Isotopes  reversed  and  the 
Isotopic  order  of  these  spectra  corresponds  to  the  Isotopic 
order  used  In  Chapter  III. 

Figure  42  shows  the  oxygen-18  concentration  versus 
depth  profile  for  a sample  oxidized  at  1473±5°K.  The  fact 
that  the  oxygen-18  "marker"  Is  located  at  the  silica- 
silicon  interface,  (Case  II),  Indicates  that  the  rate 
controlling  mechanism  for  the  thermal  oxidation  of  silicon 
at  this  temperature  Is  either  the  transport  of  molecular 
oxygen  with  the  rate  of  permeation  less  than  the  rate  of 
exchange  (see  Section  3.2.3)  or  the  transport  of  an  Ionic 
species  (see  Sections  3.3.1  and  3.3.2). 

Figure  43  shows  the  concentration  versus  depth  pro- 
file for  a sample  oxidized  at  1412±5°K.  Here  the  shape  and 
position  of  the  oxygen-18  "marker"  corresponds  to  a situation 
Intermediate  to  Cases  I and  II  (see  Section  3.2.4).  The  rate 
controlling  mechanism  at  this  temperature  is  thus  Interpreted 
to  be  the  transport  of  molecular  oxygen,  with  the  rate  of 
permeation  equal  to  or  slightly  less  than  the  rate  of  exchange. 

Figure  44  again  shows  the  oxygen-18  concentration 
versus  depth  profile  which  Is  expected  if  the  rate  controlling 
mechanism  is  the  transport  of  molecular  oxygen  with  the  ki- 
netics of  the  permeation  and  exchange  reaction  roughly  equal 


(see  Section  3.2.4).  This  sample  was  oxidized  at  1295±5*k. 


The  profile  shown  In  Figure  45  Is  for  the  sample 
oxidized  at  the  lowest  temperature*  1171±5°K.  At  this 
temperature  the  oxygen-18  "marker"  was  found  at  the  silica- 
oxygen  Interface,  (Case  I),  Indicating  the  rate  controlling 
mechanism  at  this  temperature  was  the  transport  of  molecular 
oxygen  with  the  rate  of  permeation  much  greater  than  the 
rate  of  exchange  (see  Section  3.2.2). 

The  fact  that  the  oxygen-18  "marker"  position  Is  a 
function  of  the  oxidizing  temperature  rules  out  the  possibil- 
ity that  the  rate  controlling  process  In  the  thermal  oxida- 
tion of  silicon  Is  the  transport  of  either  silicon  or  oxygen 
Ions.  If  the  reaction  were  controlled  by  the  transport  of 
an  Ionic  species  the  oxygen-18  "marker"  would  always  be 
located  at  the  silica-silicon  interface  (see  Sections  3.3.1 
and  3.3.2).  The  motion  of  the  "marker"  with  temperature 
Indicates  that  two  processes  are  Involved  In  the  thermal 
oxidation  reaction,  the  kinetics  of  each  having  different 
temperature  dependences.  These  studies  indicate  that  the 
two  processes  Involved  are  the  permeation  of  molecular  oxygen 
through  the  adherent  silica  scale  and  the  exchange  of  the 
permeating  oxygen  with  the  network  oxygen.  This  study  also 
suggests  that  the  kinetics  of  exchange  are  roughly  equal  to 
the  kinetics  of  permeation  at  1323°K  and  exceed  the  per- 
meation kinetics  at  higher  temperatures. 
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The  results  of  these  mechanistic  experiments,  when 
combined  with  the  pressure  dependence  reported  by  other 
investigators  (4,22)-  the  parabolic  rate  constant  being 
directly  proportional  to  the  oxygen  partial  pressure-  leads 
to  the  same  conclusion  arrived  at  by  Ligenza  and  Spitzer  (22), 
l.e.  the  rate  controlling  process  in  the  thermal  oxidation  of 
silicon  is  the  diffusion  of  molecular  oxygen  through  the 
adherent  silica  scale  with  appreciable  isotopic  exchange 
occurring  between  the  diffusing  oxygen  and  the  network  oxygen. 

According  to  Meek  (26)  oxygen  diffusion  in  silica 
Involves  two  processes,  as  determined  in  these  studies,  a 
transport  process  and  an  exchange  process.  If  the  mechanism 
for  the  thermal  oxidation  of  silicon  (the  transport  of 
molecular  oxygen  through  the  adherent  silica  scale)  presented 
here  is  to  be  consistant  with  Meek's  argument  it  should  be 
possible  to  calculate  the  diffusivity  of  oxygen  in  thermally 
grown  silica  using  the  parabolic  rate  constant  determined 
in  the  kinetic  studies.  This  value  should  be  approximately 
equal  to  the  diffusivity  measured  by  Norton  in  his  permeation 
experiment.  The  diffusivity  can  be  related  to  the  parabolic 
rate  constant  by  the  following  equation  (3,4): 


D 


(43) 


where  C"  Is  the  concentration  of  oxygen  incorporated  in  the 
silica  scale  (2.2  x 10  2 molecules  / cm3)  (26)  and  C'  is  the 
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solubility  of  oxygen  In  silica  (5.6  x 10  molecules  / cm  ) 

(26).  Using  these  values  one  obtains  the  following  expres- 
sion for  the  diffusivity: 

D * [1.47  x 10"4  exp  - (1-ffieV)]  f~  (44) 

This  calculated  value  Is  within  experimental  error  of  Norton's 
calculated  diffusivity  (21): 

D = [2.88  x 10‘4  exp  - t1 -1J-V)]  (45) 

Of  even  greater  interest  is  the  fact  that  the  final 
scales  In  the  successive  isotopic  oxidation  studies  should 
be  composed  of  ten  percent  oxygen-18,  while  it  is  found  that 
only  twenty-five  to  fifty  percent  of  the  oxygen-18  incor- 
porated in  the  silica  scale  during  the  initial  oxidation 
remains  at  the  end  of  the  experiment  (see  Figures  42-45). 

This  decrease  in  oxygcn-18  concentration  must  occur  by 
diffusion  of  oxygen-18  out  of  the  silica  scale  at  the  silica- 
oxygen  interface.  The  oxygen-18  cannot  be  lost  by  any  other 
diffusion  of  oxygen  in  any  other  direction  since  this 
oxygen-18  would  be  located  and  measured  during  the  proton 
activation  analysis.  Diffusion  of  oxygen  across  the  silica- 
oxygen  interface  requires  diffusion  of  oxygen  up  an  oxygen 
concentration  gradient  but  down  an  isotopic  concentration 
gradient.  Although  it  is  not  possible  to  calculate  a dif- 
fusion coefficient  from  the  "marker"  experimental  data  due 
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to  the  qualitativeness  of  the  fits.  It  Is  possible  to  make 
some  approximate  calculation.  Shewmon  (56)  gives  the  follow- 
ing equation  for  the  degassing  of  thin  metal  films  when  the 
ratio  of  the  average  final  composition  of  a sample  to  Its 
average  Initial  composition  Is  less  than  0.8: 

2 

Average  Final  Composition  _ 8 „„„  ,tir  D»  . 

Average  Initial  Composition  ~ ^ e p ” ' ' (46) 

Using  this  equation  it  is  possible  to  calculate  a range  of 
dl ffusivl ties  corresponding  to  the  amount  of  oxygen-18 
remaining  In  the  silica  after  the  second  oxidation.  The 

_ l **  _ 1 3 

calculated  diffusivlties  fall  In  the  range  10  to  10 
2 

cm  /sec  for  the  variation  in  oxygen-18  concentration 
observed.  This  range  of  diffusivlties  Is  In  agreement 
with  the  diffusivity  measured  by  Williams  (20)  In  his 
Isotopic  exchange  experiments: 

D = [2.0  x 10’9  exp  - (^eV)]  (47) 

It  should  also  be  noted  that  less  oxygen-18  remains 
in  the  silica  scale  formed  at  the  lower  temperature.  This 
can  possibly  be  attributed  to  the  fact  that  at  low  tempera- 
tures the  oxygen-18  remains  peaked  at  the  si  1 ica-oxygen 
Interface  giving  a shorter  diffusion  distance  for  the 
removal  of  oxygen-18. 

These  "marker"  studies  therefore  indicate  that  the 
diffusivity  of  oxygen  in  silica  governing  the  transport  of 
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oxygen  Is  that  value  Measured  by  Norton  (20)  In  his  per- 
meation experiments.  The  transport  of  oxygen  therefore 
occurs  by  oxygen  dissolving  In  the  silica  at  the  silica- 
oxygen  Interface  and  diffusing  through  the  oxide.  An 
exchange  also  occurs  between  the  diffusing  oxygen  and  the 
network  oxygen.  The  exchange  probably  occurs  by  dissolved 
oxygen  molecules  entering  vacancies  which  exist  in  the 
silica  scale.  At  thes*  sites  broken  slllcon-to-oxygen 
bonds  will  exist  which  makes  It  possible  for  one  of  the 
oxygens  to  become  part  of  the  oxide.  It  is  equally  likely 
that  an  oxygen  which  Is  bonded  to  only  one  silicon  can 
bond  to  the  single  oxygen  atom,  which  has  lost  an  oxygen 
atom  to  the  scale,  and  continue  diffusing  through  the 
scale.  This  model  Is  consistent  with  the  model  proposed 
by  Meek  (26)  for  oxygen  diffusion  In  vitreous  silica. 
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SECTION  VII 

SUGGESTIONS  FOR  FUTURE  WORK 


7.1  Kinetic  Oxidations 

It  Is  felt  that  since  the  thermal  oxide  scales  found  in  ‘ ] 

these  investigations  are  thicker  than  others  reported  in  the  liter- 
ature a series  of  experiments  should  be  undertaken  to  investigate 
the  effect  of  water  vapor  content  of  the  ambient  on  the  parabolic 
rate  constant.  It  is  also  felt  that  the  temperature  range  and  time 
of  oxidation  should  be  extended  in  order  to  better  determine  the 
exact  kinetics  of  the  reaction. 

7.2  Mechanistic  Oxidations 

The  mechanistic  experiment  has  shown  that  the  relative  rate 

of  exchange  becomes  approximately  equal  to  the  rate  of  permeation  around 
0 

1323UK,  and  then  becomes  more  rapid  than  the  rate  of  permeation  at  higher 
temperatures.  This  result  is  based  on  runs  at  only  four  temperatures. 

It  is  felt  that  the  number  of  temperatures  investigated  should  be  increased 
to  allow  better  determination  of  the  relative  exchange  and  permeation  rates. 

Along  the  same  line  it  is  felt  that  a numerical  solution,  and 
possibly  an  analytical  solution,  exists  to  predict  the  oxygen-18  distri- 
butions expected  in  tracer  experiments  where  the  distribution  will  be  a 
function  of  exchange  rate  and  permeation  rate. 


7.3  Experimental  Procedure 

The  construction  of  an  oxidation  apparatus  that  can  accommodate 
samples  up  to  1.5  centimeters  square  will  be  invaluable.  (This  apparatus 
is  under  construction).  The  larger  samples  will  allow  the  use  of  ellips- 
ometry  as  well  as  proton  activation  techniques. 


I 


117 


T 


In  future  studies  where  it  is  possible  to  produce  thin 
films  of  known  thickness  and  uniform  oxygen-18  concentration  the 
stopping  powers  for  protons  and  alphas  in  the  material  should  be 
checked.  This  experimental  procedure,  which  has  been  outlined  in 
Chapter  VI,  is  straight  forward  if  one  works  in  a region  of  the 
cross  section  that  is  rapidly  changing  (i.e.  the  637  KeV  resonance) 
to  give  large  changes  in  alpha  flux  for  small  changes  in  energy. 

7.4  Data  Reduction 

A deconvolution  program  for  the  analysis  of  alpha  spectra 
from  oxidation  and  diffusion  runs  would  greatly  decrease  the  amount 
of  time  required  to  analyze  each  spectrum.  The  trial  and  error 
convolution  procedure  presently  used  is  quite  laborious  whereas  a 

deconvolution  program  would  give  the  concentration  versus  depth 
profile  directly. 
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APPENDIX  A 


RANGE-ENERGY  CURVES  FOR  PROTONS  AND  ALPHAS  ON  Ta^g,  ZnFg,  AND 
Si02  AND  KINEMATICS  OF  THE  180(p,a)15N,  19F(p,a)  160,  AND 
19F(p.Ol)160*  NUCLEAR  REACTIONS 
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TABLE  A1  --Range-energy  curve  for  protons  on  Ta205 


Depth 

In 

microns 

0.00 

0.04 

0.08 

0.12 

0.16 

0.000 

1000  ocooo 

996.83937 

993.67426 

990.50311 

987.32748 

0.200 

984.14572 

98C. 95945 

977.76863 

974.57165 

971.3701C 

0.400 

968.16232 

964.94994 

961.73295 

958.50967 

955.28174 

0.600 

952.04747 

948.80845 

945.56300 

942.31287 

939.05637 

0.600 

935.79516 

932.52751 

929.25514 

925.97626 

922.69263 

1.000 

919.40246 

916.10572 

912.80417 

909.49501 

906.18300 

1.200 

902.86333 

899.53876 

896.20749 

892.86951 

889.52564 

1.400 

886.17703 

832.82064 

S79. 45744 

876.03740 

872.71239 

1.600 

869.33050 

865.94159 

862.54787 

859.14515 

855.73738 

i.eoo 

852.32260 

84S. 90067 

845.47357 

842.03754 

838.59644 

2.000 

835.14830 

831.69308 

828.23074 

824.76123 

821.28456 

2.200 

617.80066 

814.30952 

810.31110 

807.30538 

803.75232 

2.400 

800.27 189 

796.74377 

7S3. 20833 

789.66553 

786.11533 

2.600 

782.55769 

788.99040 

775.41779 

771.83543 

768.24549 

2.800 

764.64792 

761.04044 

757.42754 

753.80462 

750.17394 

3.000 

746.53526 

742.86653 

739.23002 

735.56571 

731.83116 

3.200 

728.20871 

724.51832 

780.81754 

717.10873 

713.39187 

3. 400 

709.65443 

7ac  eococ 

702.18012 

69S. 42547 

694.65999 

3.600 

690.88379 

687.09937 

682.30663 

679.5C309 

675.68352 

3.800 

671.86291 

668.02616 

664.18362 

660.32983 

656.46475 

4.000 

652.53830 

648.70021 

644.30059 

640.88958 

636.96712 

4.2C0 

633.03594 

629.09602 

625.14443 

621.18109 

617.20593 

4.400 

613.21835 

609.21979 

605.20867 

601.18542 

597.14S31 

4.600 

593.10121 

589.04097 

584.96538 

580.87429 

576.77078 

4.300 

572.65472 

558.52604 

564.38468 

560.23055 

556.06357 

5.000 

551.88026 

547.67952 

543.46641 

539.24046 

535.00155 

5.200 

530.74261 

526.47045 

522.18493 

517.88245 

513.56276 

5.400 

509.22942 

504.67858 

500.51010 

496.12646 

491.72545 

5.600 

487.30692 

482.87062 

478.41639 

473.94802 

469.45727 

5.800 

464.95211 

460.42408 

455.88138 

451.31532 

446.72891 

6.000 

442.12381 

437.49525 

432.85178 

428.18426 

423.49694 

6.200 

418.78495 

414.05274 

409.29526 

404.51708 

399.71092 

6.400 

394.88436 

390.03249 

385.15485 

380.25104 

375.32588 

6.600 

370.37382 

365.39448 

360.38748 

355.35245 

350.28905 

6.800 

345.19476 

340.06727 

334.91163 

329.72734 

324.50777 

7.000 

319.25853 

313.97231 

308.64935 

303.29575 

297.90131 

7.200 

292.46803 

287.00415 

281.49515 

275.95421 

270.36604 

7.400 

264.73716 

259.C6666 

253.34565 

247.57949 

241.76930 

7.600 

235.90397 

229.99705 

224.03216 

218.01300 

211.93828 

7.800 

2C5. 80582 

199.61748 

193.35749 

187.02050 

180.61795 

8.000 

174.13811 

167.53110 

160.93481 

154.21126 

147.41069 

8.200 

140.52064 

133.54113 

126.47221 

119.31398 

112.05337 

8.400 

104.70378 

97.25209 

89.69824 

82.02875 

74.27154 

8.600 

65.39085 

58.41067 

50.20733 

42.09004 

33.74294 

8.800 

25.26705 

16.66247 

7.92985 

.00000 

.00000 
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TABLE  A2  -Range-energy  curve  for  alphas  on  TagO,. 


Depth 

In  0.00  0.04  0.08  0.12  0.16 

Microns 


0.000 

0.200 

0.400 

0.600 

0.800 

1.000 

1.200 

1.400 
1.600 
1.800 
2.000 
2.200 

2.400 
2.600 
2.800 

3.000 

3.200 

3.400 

3.600 
3.eoo 

4.000 

4.200 

4.400 

4.600 

4.800 

5.000 

5.200 

5.400 

5.600 

5.800 

6.000 

6.200 

6.400 


5000.00000 

4943.16060 

4885.96940 

4328.42090 

4770.EC5C0 

4712.21700 

4653.55080 

4594.49570 

4535.03520 

4475.17430 

4414.89300 

4354.17970 

4293.02730 

4231.41950 

4169.34340 

4106.78820 

4043.74570 

3980.19900 

3916.13510 

3851.53440 

3786.37910 

3720.65220 

3654.33560 

3587.41070 

3519.86040 

3451.55860 

3382.78260 

3313.21000 

3242.91940 

3171.88280 

3100.06510 

3027.43420 

2953.96100 


4988.66000 
4931.75070 
4874.48900 
4816.85690 
4758.8770G 
4700.51470 
4641.77080 
4582.63610 
4523.09620 
4463.15220 
4402.73390 
4341.93590 
4280.74140 
4219.04340 
4156.87070 
4094.21930 
4031.07750 
3967.42850 
3303.25760 
3833.54810 
3773.28 020 
3707.43650 
3641.00130 
3573.95220 
3505.27240 
3437.93510 
3368.92470 
3299.20980 
3228.77370 
3157.53350 
3085.60570 
3012.60710 
2939.16160 


4977.30580 
4920.32650 
4862.99310 
4805.29840 
4747.23450 
4638.97660 
4629.97460 
4570.76010 
4511.14070 
4451.11330 
4390.65840 
4329.77210 
4263.43330 
4206.64670 
4144.37800 
4081.63090 
4018.3SS00 
3954.63650 
3390.35990 
3S25. 53960 
3760.15650 
3694.19810 
3627.64010 
3560.46650 
3492.65880 
3424.19140 
3355.03870 
3285.18080 
3214.59860 
3143.25200 
3071.11130 
2998.14670 
2924.32810 


4565.93740 
4908.88800 
4851.48390 
4793.71520 
4735.57750 
4677.06230 
4618.16470 
4558.86780 
4499.16850 
4439.05770 
4378.51580 
4317.54110 
4256.11770 
4194.23170 
4131.86790 
4069.02310 
4005.67870 
3941.82440 
3877.44020 
3612. 5uts8G 
3747.01330 
3680.93430 
3614.25560 
3546.95720 
3479.01930 
3410.41540 
3341.12430 
3271.12280 
3200.39170 
3128.88320 
3056.58620 
2983.45340 
2909.45790 


4954.55610 

4897.43650 

4839.9590C 

4782.11730 

4723.90460 

4665.31460 

4606.33870 

4546.95900 

4487.17970 

4425.93520 

4366.33630 

4335.29290 

4243.77810 

4181.79700 

4119.33760 

4056.39420 

3592.94960 

3928.99070 

3864.45840 

3799.45560 

3733.84500 

3667.64680 

3600.84580 

3533.42060 

3465.35200 

3395.61240 

3327.18140 

3257.03570 

3186.15210 

3114.49400 

3042.02560 

2963.72410 

2894.55030 


TABLE  A2  —Continued 


Depth 

In 

microns 

0.00 

0.04 

0.08 

0.12 

0.16 

6.600 

2879.50510 

2864.62310 

2849.60410 

2834.54800 

2819.45450 

6.800 

2804.32360 

2789.15440 

2773.94670 

2758.70110 

2743.41470 

7.000 

2728.03720 

2712.72130 

2697.31410 

2681.85420 

2666.37340 

7.200 

2650.84150 

2635.26840 

2619.65370 

2603.99440 

2588.29640 

7.400 

2572.55300 

2556.76400 

2540.92920 

2525.05130 

2509.130C0 

7.600 

2493.16210 

2477.14730 

2461.08550 

2444.97660 

2428.81450 

7.800 

2412.60480 

2396.34740 

2380.04240 

2363.68320 

2347.275S0 

8.0C0 

2330.81660 

2314.30560 

2297.74250 

2281.12650 

2264.45490 

8.200 

2247.73070 

2230.95360 

2214.12020 

2197.23000 

2180.28310 

8.400 

2163.28220 

2146.22350 

2129.10340 

2111.92510 

2094.68820 

8.600 

2077.38360 

2060.02990 

2042.60440 

2025.11930 

2007.56710 

8.800 

1989.95530 

1972.27670 

1954.53100 

1936.72510 

1918.85130 

9.000 

1900.90910 

1882.89610 

1864.81510 

1846.66520 

1828.44620 

9.200 

1810.14960 

1791.78330 

1773.34240 

1754.82540 

1736.23200 

9.400 

1717.56180 

1698.81020 

1679.98100 

1661.06990 

1642.07660 

9.600 

1623.0C070 

1603.84170 

1584.59930 

1565.26340 

1545.84880 

9.800 

1526.33590 

1506.74150 

1457.04870 

1467.27040 

1447 . 33G6C 

10.000 

1427.42700 

1407.36130 

1387.19940 

1366.94060 

1346.57470 

10.200 

1326.11160 

1305.54100 

1284.87410 

1264.09840 

1243.21420 

10.40U 

1222.22160 

1201.12070 

1179.91190 

1158.58910 

1137.1579C 

10.600 

1115.61340 

1093.95610 

1072.17970 

105C. 28430 

1028.27080 

10.800 

1006.14029 

983.89447 

961.53241 

939.04511 

916.44385 

11.000 

893.72102 

870.88324 

847.92724 

824.85508 

801.67031 

11.200 

778.37332 

754.96515 

731.45659 

707.84472 

654.14058 

11.400 

660.34794 

636.47360 

612.52424 

588.51250 

564.44987 

11.600 

540.35441 

516.23589 

492.12109 

468.02177 

443.97245 

11.800 

419.99701 

396.12010 

372.39090 

348.84964 

325.53982 

12.000 

302.53191 

279.85414 

257.56913 

235.78061 

214.55834 

12.200 

194.00313 

174.11718 

155.06183 

136.95589 

119.92844 

12.400 

104.08955 

89.43985 

76.03236 

64.01065 

53.06524 

12.600 

43.33601 

34.60065 

26.92987 

20.10485 

14.13502 

12.800 

8.88307 

4.32130 

.26947 

.00000 

.00000 

t 

1 
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TABLE  A3  —Kinematics  of  the  180(p,a)15N  reaction 


PROTON 

ALPHA 

PROTON 

ALPHA 

ENERGY 

ENERGY 

ENERGY 

ENERGY 

(MeV) 

(MeV) 

(MeV) 

(MeV) 

0.500 

3.250 

0.760 

3.374 

0.510 

3.254 

0.770 

3.379 

0.520 

3.259 

0.780 

3.384 

0.530 

3.263 

0.790 

3.389 

0.540 

3.268 

0.800 

3.394 

0.550 

3.273 

0.810 

3.399 

0.560 

3.277 

0.820 

3.404 

0.570 

3.282 

0.830 

3.409 

0.580 

3.287 

0.840 

3.414 

0.590 

3.292 

0.850 

3.419 

0.600 

3.296 

0.860 

3.424 

0.610 

3.301 

0.870 

3.430 

0.620 

3.306 

0.880 

3.435 

0.630 

3.310 

0.890 

3.440 

0.640 

3.315 

0.900 

3.445 

0.650 

3.320 

0.910 

3.450 

0.660 

3.325 

0.920 

3.455 

0.670 

3.330 

0.930 

3.460 

0.680 

3.335 

0.940 

3.465 

0.690 

3.340 

0.950 

3.470 

0.700 

3.345 

0.960 

3.476 

0.710 

3.350 

0.970 

3.481 

0.720 

3.354 

0.980 

3.486 

0.730 

3.359 

0.990 

3.491 

0.740 

3.364 

1.000 

3.496 

0.750 

3.369 
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TABLE  A4  -Range-energy  curve  for  protons  on  ZnF2 


Depth 

In  0.00  0.04  0.08  0.12  0.16 

■Icrons 


0.000 

1000.00000 

977.19401 

994.38459 

991.57020 

988.75079 

0.200 

985.92792 

983.10003 

980.26708 

977.43066 

974.58915 

0.400 

971.74256 

968.89248 

966.03728 

963.17693 

960.31307 

0.600 

957.44406 

954.55985 

951.69047 

948.80591 

945.91781 

0.800 

943.02452 

940.12599 

937.22221 

934.31483 

931.40230 

1.  x»o 

928.48441 

925.56123 

922.63275 

919.69892 

916.75974 

1.400 

913.81696 

910.86879 

907.91524 

904.95629 

901.99192 

1.400 

899.02199 

896.04652 

893.06565 

890.07934 

887.08761 

1.600 

884.09042 

881.08775 

878.07961 

875.06595 

872.04678 

1.800 

869.02206 

865.99178 

862.95594 

859.91450 

856.86745 

2.000 

853.81479 

850.75648 

847.69056 

844.61901 

841.54180 

2.200 

838.45892 

835.37035 

832.27608 

829.17410 

826.06638 

2.400 

822.95290 

819.83365 

816.70860 

813.57572 

810.43700 

2.600 

807.29244 

804.13996 

800.98160 

797.81704 

794.64668 

2.800 

791.46840 

788.28423 

735.09205 

781.89395 

778.68991 

3.000 

775.47776 

772.25964 

769.03334 

765.80102 

762.56047 

3.200 

759.31335 

756.05893 

752.79791 

749.52832 

746.25259 

3.400 

742.96856 

739.67844 

736.37996 

733.07305 

729.75998 

3.600 

726.43842 

723.10831 

719.77196 

716.42699 

713.07337 

3.800 

709.71344 

706.34478 

702.96735 

699.58322 

696.19025 

4.000 

692.78848 

689.37791 

685.95843 

682.53263 

679.09533 

4.200 

675.65156 

672.19875 

668.73687 

665.26587 

661.78571 

4.400 

658.29634 

654.79772 

551.28981 

647.77233 

644.24556 

4.600 

640.70969 

637.16438 

633.60965 

630.04274 

626.4690? 

4.800 

622.88311 

619.28754 

615.68233 

612.06456 

608.43994 

5.000 

604.80265 

601.15555 

597.49794 

593.82775 

590.14779 

5.200 

586.45796 

582.75524 

579.03953 

575.31685 

571.58104 

5.400 

567.83200 

564.07279 

560.30336 

556.52048 

552.72406 

5.600 

548.91384 

545.09316 

541.26219 

537.41764 

533.55938 

5.800 

529.68733 

525.80138 

521.90144 

517.98739 

514.06612 

6.000 

510.13061 

506.18076 

502.21648 

498.23718 

494.24312 

6.200 

490.23481 

486.21213 

482.17496 

478.12318 

474.05667 

6.400 

469.97148 

465.86746 

461.74830 

457.61390 

453.46414 

1'  6.600 

449.29804 

445.11632 

440.91532 

436.69487 

432.45883 

6.800 

428.20710 

423.93953 

419.64740 

415.33915 

411.01465 

7.000 

406.66483 

402.29847 

397.91494 

393.50547 

389.07957 

7.200 

384.62777 

380.15918 

375.66885 

371.15660 

366.62214 

7.400 

362.06531 

357.48072 

352.87847 

348.24689 

343.59183 

1 7.600 

338.91430 

334.20867 

329.47985 

324.72209 

319.94060 

7.800 

315.12938 

310.28801 

305.42199 

300.52508 

295.59452 

8.000 

290.63407 

285.64288 

280.62052 

275.56645 

270.48019 

8.200 

265.36123 

260.20913 

255.01646 

249.78966 

244.51876 

1 8.400 

239.21416 

233.86799 

228.47933 

223.05524 

217.57936 

| 8.600 

212.06696 

206.50098 

200.89751 

195.24211 

189.52146 

6.800 

183.75685 

177.92827 

172.03648 

166.08223 

160.06625 

. 9.000 

9.200 

153.98926 

147.85197 

141.65510 

135.38917 

129.05505 

122.66362 

116.20543 

109.68132 

103.09202 

96.43833 

9.400 

89.71096 

82.92077 

76.05852 

69.12507 

62.12127 

9.600 

55.04792 

47.90586 

40.69589 

33.40895 

26.04591 

9.800 

18.61753 

11.10493 

3.52873 

.00000 

.00000 
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TABLE  A5  -Range-energy  curve  for  alphas  on  ZnF2 


Death 


In 

•1 crons 

0.00 

0.04 

0.08 

0.12 

0.16 

0.000 

5000.00000 

4990.16020 

4980.30860 

4970.44510 

4960.57100 

0.200 

4950.68620 

4940.73950 

4930.88100 

4920.96060 

4911.02810 

0.400 

4901.08380 

4891.12740 

4881.15910 

4871.17860 

4861.18620 

0.600 

4851.18170 

4841.16510 

4831.13640 

4821.09560 

4811.04260 

0.800 

4800.97740 

4790.89980 

4780.81000 

4770.70810 

4760.59390 

1.000 

4750.46630 

4740.32520 

4730.17190 

4720.00610 

4709.82820 

1 

1.200 

4699.63790 

4689.43390 

4679.21630 

4668.98630 

4658.74380 

! 

1.400 

4648.43880 

4638.21890 

4627.93650 

4617.64150 

4607.33260 

1.600 

4597 . 0098C 

4536.67430 

4576.32500 

4565.96170 

4555.58600 

1.800 

4545.19500 

.4534.79130 

4524.37370 

4513.94210 

4503.49780 

2.000 

4493.03810 

4482.56550 

4472.07750 

4461.57670 

4451.06C40 

2.200 

4440.53110 

4429.98630 

4419.42860 

4408.85510 

4398.26710 

2.400 

4387.66460 

4377.04780 

4366.41660 

4355.76970 

4345.10830 

2.600 

4334.43240 

4323.74060 

4313.03430 

4302.31200 

4291.57500 

2.800 

4280.82210 

4270.05450 

4259.27080 

4248.47250 

4237.65790 

3.000 

4226.82700 

4215.98140 

4205.11960 

4194.24000 

4183.34530 

3.200 

4172.43580 

4161.51150 

4150.57210 

4139.61760 

4128.64780 

3.400 

4117.66260 

4106.66180 

4095.54510 

4084.61250 

4073.56400 

3.600 

4062.49920 

4051.41610 

4040.32050 

4029.20480 

4018.07230 

3.800 

4006.92300 

3995.75580 

3984.56690 

3973.35820 

3962.12790 

1 

4.000 

3950.87790 

3939.60820 

3928.31720 

3917.00480 

3905.67310 

i 

4.200 

3894.32020 

3882.94610 

3871.55290 

3860.13690 

3848.70190 

4.400 

3237.24610 

3325.76340 

3814.27220 

3S02. 75440 

3791.23630 

4.600 

3775.69560 

3763.13540 

3756.55260 

3744.94630 

3733.32170 

4.800 

3721.67330 

3710.00190 

3698.31120 

3686.59770 

3674.86140 

i 

5.000 

3663.10420 

3651.32610 

3639.52550 

3627.70240 

3615.85680 

5.200 

3603.98880 

3592.08900 

3580.15760 

3568.20880 

3556.24240 

5.400 

3544.25830 

3532.25610 

3520.23570 

3508.19700 

3496.13960 

5.600 

3484.06340 

3471.96670 

3459.84910 

3447.71220 

3435.55560 

5.800 

3423.37920 

3411.18280 

3398.96630 

3386.73450 

3374.48010 

6.000 

3362.20290 

3349.90130 

3337.57480 

3325.22550 

3312.85130 

6.200 

3300.45220 

3288.03000 

3275.58070 

3263.10830 

3250.60860 

6.400 

3238.08570 

3225.53520 

3212.96140 

3200.36010 

3167.73280 

6.600 

3175.08010 

3162.39980 

3149.69380 

3136.96210 

3124.20260 

!l 

6.800 

3111.41710 

3098.60370 

3085.76200 

3072.89410 

3059.99790 

7.000 

3047.07330 

3034.12010 

3021.13830 

3008.12780 

2995.08840 

J 

7.200 

2982.02010 

2968.92290 

2955.79650 

2942.63880 

2929.45190 

m 


TABLE  A5  Continued 


! 


Depth 

In 

*1  crons 

0.00 

0.04 

0.08 

0.12 

0.16 

7.400 

2916.23560 

2902.98770 

2889.71040 

2876.40120 

2863.06250 

7.600 

2849.69180 

2336.28390 

2822.85380 

2309.38630 

27S5. 86850 

7.800 

2782.35580 

2768.79060 

2755  19300 

2741.56010 

2727.89690 

8.000 

2714.19830 

2700.46400 

2686.69880 

2672.S5770 

2659.05040 

8.200 

2645.18690 

2631.27710 

2617.33060 

2603.34770 

2589.32840 

8.400 

2575.27190 

2561.17850 

2547.04770 

2532.87960 

2518.67410 

8.600 

2504.42830 

2490. 14 ISO 

2475.81760 

2451.45570 

2447.05040 

8.800 

2432.60700 

2418.11990 

2403.59460 

2389.02540 

2374.41730 

9.000 

2359.76580 

2345.07230 

2330.33720 

2315.55750 

2300.73290 

9.200 

2285.86590 

2270.95430 

2255.99790 

2240.99630 

2225.94940 

9.400 

2210.85700 

2195.71860 

2180.53110 

2165.2SS1C 

2150.01610 

9.600 

2134.68780 

2119.31010 

2103.88240 

2038.40790 

2072.87890 

9.800 

2057.29840 

2041.66540 

2025.97930 

2010.24390 

1994.45400 

10.000 

1978.60960 

1962.71010 

1946.75510 

1930.74460 

1914.57320 

10.200 

1898.55570 

1882.37450 

1866.13780 

1349.34150 

1833.48570 

10.400 

1817.07320 

1800.5967C 

1784.06510 

1767.46730 

1750.80690 

10.600 

1734.08390 

1717.29450 

1700.44220 

1683.52310 

1665.53310 

10.800 

1649.48720 

1632.36990 

1615.18600 

1597.93110 

1580.60380 

11.000 

1563.20910 

1545.74280 

1523.20450 

15 iO. 59410 

1492.91230 

11.200 

1475.15360 

1457.31660 

1439.40110 

1421.40710 

1403.33480 

11.400 

1385.18420 

1356.95510 

1348.64750 

1330.25720 

1311.78410 

11.600 

1293.22860 

1274.59040 

1255.36490 

1237.05660 

1213.16140 

11.800 

1199.17920 

118C. 11060 

1160.95450 

1141.71150 

1122.37730 

12.000 

1102.95240 

1083.44280 

1063.83720 

1044.14130 

1024. 35500 

12.200 

1004.47764 

984.50726 

964.44520 

944.29353 

924.05310 

12.400 

903.71614 

883.29293 

862.7802 8 

842.18204 

821.51028 

12.600 

800.73308 

779.88632 

758.96055 

737.96099 

716.88791 

12.800 

695.74608 

674.54001 

653.27863 

631.96689 

610.61110 

13.000 

589.22438 

567.80986 

546.37941 

524.94670 

503.52357 

13.200 

482.13010 

460.77687 

439.49234 

418.28785 

397.19035 

13.400 

376.23817 

355.44199 

334.86611 

314.50749 

294.43271 

13.600 

274.66574 

255.24753 

236.27876 

217.73515 

199.69492 

13.800 

182.27444 

165.42606 

149.27131 

133.93268 

119.40576 

14.000 

105.84727 

93.21437 

81.50187 

70.77483 

60.91091 

14.200 

51.94603 

43.91575 

36.52927 

29.99504 

24.09231 

14.400 

18.79689 

14.07300 

9.88054 

6.17252 

2.88440 
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TABLE  A6 


19  16 

—Kinematics  of  the  F(p,a)  0 reaction 


PROTON 

ENERGY 

(MeV) 

ALPHA 

ENERGY 

(MeV) 

PROTON 

ENERGY 

(MeV) 

ALPHA 

ENERGY 

(MeV) 

0.500 

6.516 

0.760 

6.625 

0.510 

6.519 

0.770 

6.630 

0.520 

6.523 

0.780 

6.634 

0.530 

6.527 

0.790 

6.639 

0.540 

6.531 

0.800 

6.643 

0.550 

6.535 

0.810 

6.648 

0.560 

6.539 

0.820 

6.652 

0.570 

6.544 

0.830 

6.657 

0.580 

6.548 

0.840 

6.662 

0.590 

6.552 

0.850 

6.666 

0.600 

6.556 

0.860 

6.671 

0.610 

6.560 

0.870 

6.675 

0.620 

6.564 

0.880 

6.680 

0.630 

6.568 

0.890 

6.685 

0.640 

6.573 

0.900 

6.689 

0.650 

6.577 

0.910 

6.694 

0.660 

6.581 

0.920 

6.699 

0.670 

6.586 

0.930 

6.704 

0.680 

6.590 

0.940 

6.708 

0.690 

6.594 

0.950 

6.713 

0.700 

6.599 

0.960 

6.718 

0.710 

6.603 

0.970 

6.723 

0.720 

6.607 

0.980 

6.727 

0.730 

6.612 

0.990 

6.732 

0.740 

6.616 

1.000 

6.737 

0.750 

6.621 

TABLE  A7  -Kinematics  of  the  l9F(ptOj)160*  reaction 


PROTON 

ENERGY 

(MeV) 

ALPHA 

ENERGY 

(MeV) 

PROTON 

ENERGY 

(M?V) 

ALPHA 

ENERGY 

(MeV) 

0.500 

1.773 

0.760 

1.917 

0.510 

1.779 

0.770 

1.922 

0.520 

1.784 

0.780 

1.928 

0.530 

1.790 

0.790 

1.934 

0.540 

1.795 

0.800 

1.939 

0.550 

1.801 

0.810 

1.945 

0.560 

1.806 

0.820 

1.951 

0.570 

1.811 

0.830 

1.956 

0.580 

1.817 

0.840 

1.962 

0.590 

1.822 

0.850 

1.968 

0.600 

1.828 

0.860 

1.973 

0.610 

1.833 

0.870 

1.979 

0.620 

1.839 

0.880 

1.985 

0.630 

1.844 

0.890 

1.990 

0.640 

1.850 

0.900 

1.996 

0.650 

1.856 

0.910 

2.002 

0.660 

1.861 

0.920 

2.007 

0.670 

1.867 

0.930 

2.013 

0.680 

1.872 

0.940 

2.019 

0.690 

1.878 

0.950 

2.024 

0.700 

1.883 

0.960 

2.030 

0.710 

1.889 

0.970 

2.036 

0.720 

1.894 

0.980 

2.041 

0.730 

1.900 

0.990 

2.047 

0.740 

1.906 

1.000 

2.053 

0.750 

1.911 
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TABLE  A8  -Range-energy  curve  for  protons  on  Si02 


Depth 

In 

ailerons 

0.00 

0.04 

0.08 

0.12 

0.16 

0.000 

1000.00000 

998.27081 

996.53918 

994.80727 

993.07290 

0.200 

991.33608 

989.59787 

987.85828 

986.11622 

984.37169 

0.400 

982.62577 

980.87846 

979.12665 

977.37634 

975.62263 

0.600 

973.86753 

972.10991 

970.34975 

968.58820 

966.82523 

0.800 

965.05973 

963.29163 

961.52107 

959.75019 

957.97675 

1.000 

956.20074 

954.42215 

952.64214 

950.86071 

949.07640 

1.200 

947.28957 

945.50016 

943.71050 

941.91824 

940.12338 

1.400 

938.32592 

S36. 52585 

934.72554 

932.92261 

931.11705 

1.600 

929.303S5 

927.49301 

925.68693 

923.87321 

922.05682 

1.800 

920.23778 

913.41606 

916.59290 

914.76828 

912.94098 

2.000 

911.11037 

SOS. 27827 

907.44236 

905.60599 

9C3. 76767 

2.200 

9C1. 92662 

900.08291 

398.23649 

896.38734 

894.53544 

2.400 

892.63336 

890.82353 

888.97094 

887.11060 

885.24748 

2.600 

883.38160 

881.51293 

879.64408 

877.77244 

875.89801 

2.  SCO 

874.02079 

872.14074 

870.25788 

858.37219 

866.48368 

3.000 

864.59502 

862.70352 

860.80917 

858.91198 

857.01192 

3.200 

855.10901 

853.20324 

851.29458 

849.38266 

847.46802 

3.400 

845.55189 

843.63428 

841.71379 

839.79043 

837.86420 

3.600 

835.93506 

834.00305 

832.06813 

830.13030 

828.18955 

3.800 

826.24588 

824.29923 

822.34974 

820.39726 

818.44182 

4.000 

816.48342 

814.52351 

812.56209 

810.59770 

808.63033 

4.200 

806.65993 

804.68664 

802.71029 

800.73095 

798.74854 

4.400 

796.76315 

794.7747C 

792.78340 

79C. 78902 

788.79362 

4.600 

786.79120 

784.78775 

782.78127 

780.77172 

778.75914 

4.800 

776.74350 

774.72478 

772.70300 

770.67314 

768.65018 

5.000 

766.61913 

764.58498 

762.54772 

760.50735 

758.46066 

5.200 

756.41083 

754.35786 

752.30175 

750.24212 

748.17936 

5.400 

746.11349 

744.04451 

741.97239 

739.89716 

737.81379 

5.600 

735.73726 

733.65257 

731.56473 

729.47203 

727.37449 

5.800 

725.27374 

723.16979 

721.06262 

718.95223 

716.83861 

6.000 

714.72176 

712.60166 

710.47831 

708.34321 

706.21435 

6.200 

704.07819 

701.93826 

699.79487 

697.64816 

695.49819 

6.400 

693.34140 

691.18132 

689.01794 

686.85126 

684.68127 

6.600 

682.50795 

680.32765 

678.14399 

675.95699 

673.76661 

6.800 

671.57287 

669.37388 

667.16962 

664.96195 

662.75088 

7.000 

660.53640 

658.31465 

656.08947 

553.86084 

651.62876 

7.200 

649.39112 

647.14794 

644.90134 

642.65130 

640.39584 

7.400 

638.13494 

635.87058 

633.60274 

631.32939 

629.05053 

7.600 

626.76816 

624.48227 

622.18875 

619.89170 

617.59110 
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TABLE  A8-*Continued 


Oepth 


In 

■Icrons 

0.00 

0.04 

0.08 

0.12 

0.16 

7.800 

615.28486 

612.97296 

61u. 65749 

608.33631 

606.00943 

8.000 

603.67895 

601.34270 

599.00038 

596.65459 

594.30306 

8.200 

591.9*576 

589.58487 

587.21597 

584.84345 

582.46730 

8.400 

580.08303 

577.69512 

575.30129 

572.90152 

570.49809 

8.600 

568.08636 

565.67094 

563.24950 

560.82201 

558.3S845 

8.800 

555.94881 

553.50544 

551.05355 

548.59722 

546.13269 

9.000 

543.66450 

541.19027 

538.70597 

536.22355 

533.73100 

9.200 

531.23229 

528.72739 

526.21627 

523.69890 

521.17524 

9.400 

518.64530 

516.10902 

513.56640 

511.01479 

508.45939 

9.600 

505.89491 

503.32397 

500.74655 

498.16175 

495.57066 

9.800 

492.97056 

490.36670 

487.75373 

485.13427 

482.50829 

10.000 

479.87302 

477.23115 

474.58268 

471.92473 

469.26009 

10.200 

466.58874 

463.90775 

461.21997 

458.52540 

455.82103 

10.400 

453.10978 

450.39165 

447.66297 

444.92454 

442.17926 

10.600 

439.42720 

436.66516 

433.89309 

431.11406 

428.32803 

10.800 

425.53179 

422.72528 

419.91165 

417.08761 

414.25309 

11.000 

411.41139 

408.55908 

405.69613 

402.82245 

399.94047 

11.200 

397.04774 

394.14447 

391.23056 

388.30942 

385.37750 

11.400 

382.43471 

379.48098 

376.51623 

373.54038 

370.55334 

11.600 

367.55506 

364.54546 

361.52446 

358.49199 

355.4*799 

11.800 

352.39238 

349.32429 

346.24425 

343.14895 

340.04217 

12.000 

336.92383 

333.79386 

330.6521* 

327.49*47 

32*. 32488 

12.200 

321.13899 

317. 94 ICO 

314.73082 

311.50398 

308.26478 

12.400 

305.00861 

301.73994 

298.45097 

295.14657 

291.82998 

12.600 

288.49651 

285.14592 

281.78278 

278.40223 

275.00404 

12.800 

271.59299 

268.16402 

264.71692 

261.25155 

257.76771 

13.000 

254.26527 

250.74404 

247.20229 

243.64193 

240.05730 

13.200 

236.45927 

232.83648 

229.19440 

225.53282 

221.84572 

13.400 

218.13884 

214.41200 

210.65893 

206.87935 

203.07941 

13.600 

199.26359 

195.41913 

191.54222 

187.63327 

183.69263 

13.800 

179.72068 

175.71778 

171.68428 

167.62053 

163.52689 

14.000 

159.40367 

155.25123 

151.06991 

146.85253 

142.59956 

14.200 

138.31881 

134.01061 

129.66055 

125.28381 

120.87338 

14.400 

116.42967 

111.96028 

107.45113 

102.91702 

98.34397 

14.600 

93.73956 

89.10415 

84.43109 

79.72782 

74.99470 

14.800 

70.22516 

65.42652 

60.59228 

55.72287 

50.82549 

15.000 

45.89374 

40.92802 

35.92875 

30.89635 

25.83122 

15.200 

20.73374 

15.59781 

10.43040 

5.32187 

.00000 
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TABLE  A9  -Range-energy  curve  for  alphas  on  Si 0^ 


Depth 

In 

microns 


0.000 

0.200 

0.400 

0.600 

0.800 

1.000 

1.200 

1.400 
1.600 
1.800 
2.000 
2.200 

2.400 
2.600 
2.800 

3.000 

3.200 

3.400 

3.600 

3.800 

4.000 

4.200 

4.400 

4.600 

4.800 

5.000 

5.200 

5.400 

5.600 

5.800 

6.000 

6.200 

6.400 

6.600 

6.800 
7.000 
7.200 

7.400 
7.600 


5000.00000 

4970.19240 

4940.26590 

4910.22010 

4580.05400 

4849.76700 

4819.35610 

4788.81620 

4758.14860 

4727.35250 

4696.42530 

4665.36540 

4634.17010 

4602.83860 

4571.36850 

4539.75850 

4508.00550 

4476.10800 

4444.06300 

4411.86880 

4379.52170 

4347.01940 

4314.36230 

4281.54740 

4248.57140 

4215.43110 

4182.12420 

4148.64920 

4115.00110 

4081.17850 

4047.17790 

4012.99600 

3978.62800 

3944.07380 

3909.32910 

3874.39210 

3839.25810 

3803.91960 

3768.37610 


4994.04800 
4954.21660 
4934.26640 
4904.19650 
4374.00630 
4843.69520 
4813.25810 
4782.69330 
4751.99950 
4721.17710 
4690.22350 
4659.13690 
4627.91470 
4596.55620 
4565.05750 
4533.42020 
4501.63800 
4469.711C0 
4437.63640 
4405.41240 
4373.03310 
4340.50100 
4307.81190 
4274.96460 
4241.95600 
4208.78310 
4175.44310 
4141.93290 
4108.25040 
4074.39310 
4040.35670 
4006.13740 
3971.73170 
3937.14000 
3902.35740 
3867.38110 
3832 . 20640 
3796.82760 
3761.24290 


4988.09110 
495S. 23610 
4920.26190 
4898.16810 
4857.95380 
4837.61830 
4807.15520 
4776.56550 
4745.34560 
4714.99680 
4684.01670 
4652.90340 
4621.65440 
4590.26710 
4558.74120 
4517.07490 
4495.26430 
4463.30790 
4431.20380 
4398.94850 
4366.53940 
4333.97540 
4301.25510 
4268.37650 
4235.33430 
4202.12860 
4168.75440 
4135.21010 
4101.49320 
4067.60010 
4033.52770 
3999.27080 
3964.82880 
3930.19840 
3895.37790 
3860.36220 
3825.14660 
3789.72760 
3754.10050 


4982.12960 

4952.25080 

4922.25250 

4892.13490 

4861.89640 

4831.53670 

4801.04740 

4770.43200 

4739.68670 

4708.81160 

4677.80490 

4646.66490 

4615.38810 

4583.97300 

4552.41990 

4520.72450 

4488.88460 

4456.89870 

4424.76490 

4392.47930 

4360.03860 

4327.44460 

4294.69210 

4261.78110 

4228.70620 

4195.46750 

4162.05940 

4128.48100 

4094.72850 

4060.79960 

4026.69110 

3992.39750 

3957.91830 

3923.24900 

3888.39050 

3853.33550 

3818.07890 

3782.61840 

3746.95010 


4976.16340 

4946.26080 

4916.23890 

4886.09690 

4855.83420 

4825.44930 

4794.93420 

4764.29280 

4733.52290 

4702.62140 

4671.58810 

4640.42050 

4609.11590 

4577.67370 

4546.09180 

4514.36800 

4482.49970 

4450.48440 

4418.31990 

4386.00310 

4353.53260 

4320.90650 

4286.12300 

4255.17940 

4222.07180 

4188.79900 

4155.35810 

4121.74430 

4087.95720 

4053.99260 

4019.94680 

3985.51660 

3950.99990 

3916.29300 

3881.39520 

3846.30080 

3811.00320 

3775.50130 

3739.79170 


131 


Depth 

• 

in 

microns 

0.00 

0.04 

0.08 

0.12 

0.16 

7.800 

3732.62510 

3725.44930 

3718.26530 

3711.07210 

3703.87070 

8.000 

3696.66120 

3689.44240 

3682.21540 

3674.97690 

3667.73430 

8.200 

3660.43020 

3653.21790 

3645.94620 

3638.66610 

3631.37660 

8.400 

3624.07770 

3616.77040 

3609.45350 

3602.12710 

3594.79190 

8.600 

3587.44730 

3580.09310 

3572.7306C 

3565.35650 

3557.97680 

8.800 

3550.53660 

3543.18710 

3535.77760 

3528.35380 

3520.93020 

9.000 

3513.49310 

3506.04500 

3498.56840 

3491.12080 

3483.64460 

9.200 

3476.15880 

3468.66300 

3461.15750 

3453.64210 

3446.11550 

9.400 

3438.58030 

3431.03390 

3423.47890 

3415.91260 

3408.33780 

9.600 

3400.75150 

3393.15490 

3385.54830 

3377.93050 

3370.30410 

9.600 

3362.66640 

3355.01750 

3347.35980 

3339.69090 

3332.01050 

10.000 

3324.32140 

3316.62080 

3308.90880 

3301.18660 

3293.45430 

10.200 

3285.71050 

3277.95510 

3270.19080 

3262.41490 

3254.62740 

10.400 

3246.82820 

3239.01710 

3231.19730 

3223.36570 

3215.52230 

10.600 

3207.66720 

3199.80010 

3191.92060 

3184.03070 

3176.13030 

10.800 

3168.21800 

3163.29380 

3152.35760 

3144.40960 

3136.44950 

11.000 

3128.47730 

3120.45310 

3112.49690 

3104.48860 

3096.46810 

11.200 

3038.43550 

3080.39070 

3072.33380 

3064.26470 

3056.18320 

11.400 

3048.08960 

3039.98370 

3031.96240 

3023.72870 

3015.58270 

11.600 

3007.424*0 

2999.25360 

2931.07100 

2982.87440 

2574.66370 

11.800 

2966.44050 

2958.20470 

2949.95660 

2941.69250 

2933.41600 

12.000 

2925.12650 

2515.62340 

2003.50560 

2900.17520 

2391.82950 

12.200 

26S3. 46950 

2675.09630 

2b6t>.  70520 

2656.30330 

264;.  832*0 

12.400 

2841.46200 

2833.01630 

2624.55940 

2816.08720 

2807.59880 

12.600 

2799.09740 

2790.58060 

2782.05090 

2773.50460 

2764.94530 

12.800 

2755.36950 

2747.77880 

2739.17330 

2730.55290 

2721.91760 

13.000 

2713.26560 

2704.59850 

2695.91580 

2687.21600 

2678.50140 

13.200 

2669.77010 

2661.02400 

2652.26110 

2643.46320 

26 34. 6f 860 

13.400 

2625.87890 

2617.05230 

2608.20870 

2599.34810 

2590.47270 

13.600 

2581.58020 

2572.67040 

2563.74330 

2554.80070 

2545.64080 

1 

13.80u 

2536.86350 

2527.86870 

2518.85640 

2509.82640 

25OO.778S0 

14.000 

2491.71120 

2482.62590 

2473.52290 

2464.40220 

2455.26370 

14.200 

2446.10550 

2436.92930 

2427.73520 

2418.52110 

2409.28910 

14.400 

2400.03700 

2390.76720 

2381.47710 

2372.16890 

2362.84030 

14.600 

2353.49130 

2344.12410 

2334.73630 

2325.32800 

2315.89910 

14.800 

2306.44960 

2296.98120 

2287.43980 

2277.97990 

2268.44690 

15.000 

2258.89530 

2249.32070 

2239.72510 

2230.10660 

2220.46880 

15.200 

2210.81020 

2201.12820 

2191.42320 

2181.69700 

2171.94930 

15.400 

2162.17800 

2152.38300 

2142.56660 

2132.72870 

2122.86700 
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TABLE  A9  -Continued 


Depth 


In 

■Icrons 

0.00 

0.04 

0.08 

0.12 

0.16 

15.600 

2112.98120 

2103.07170 

2093.13750 

2083.17950 

2073.19750 

15.800 

2053.19140 

2053.16100 

2043.10660 

2033.02790 

2022.92490 

16.000 

2012.79770 

2002.64610 

1992.47010 

1982.26500 

1972.03560 

16.200 

1961.78180 

1951.50370 

1941.19590 

1930.8636C 

1920.50670 

16.400 

1910.11990 

1899.70830 

1889.26830 

1878.80110 

1868.30650 

16.600 

1857.78450 

1847.23500 

1836.65770 

ie26. 05280 

1815.42010 

16.800 

1804.75670 

1794.06560 

1783.34650 

1772.59620 

1761.81470 

17.000 

1751.00450 

1740.16290 

1729.29260 

1718.39070 

1707.45720 

17.200 

1696.49180 

1685.49460 

1674.46550 

1663.40440 

1652.31130 

17.400 

1641.18320 

1630.02280 

1618.83030 

1607.60540 

1596.34540 

17.600 

1585.04980 

1573.72150 

1562.35750 

1550.95760 

1539.52490 

17.800 

1528.05330 

1516.54880 

1505.00530 

1493.42900 

1481.81380 

18.000 

1470.16240 

1458.47480 

1446.74790 

1434.98170 

1423.18230 

18.200 

1411.34360 

1399.46530 

1387.54700 

1375.58900 

1363.59130 

18.400 

1351.55380 

1339.47660 

1327.35950 

1315.20270 

1303.00620 

18.600 

1290.76430 

1278.48210 

1266.15970 

1253.79380 

1241.38440 

18.800 

1228.93500 

1216.43890 

1203.90290 

1191.32030 

1178.69470 

19.000 

1166.02610 

1153.31110 

1140.55320 

1127.74920 

1114.90230 

19.200 

1102.00950 

1089.07030 

1076.08430 

1063.05190 

1049.97320 

19.400 

1036.84860 

1023.67486 

1010.45576 

997.19271 

983.88069 

19.600 

970.51964 

957.11373 

943.65977 

930.15846 

916.60731 

19.800 

903.01331 

889.37240 

875.68635 

861.95632 

848.17566 

20.000 

834.35620 

820.48672 

806.57870 

792.62781 

778.63380 

20.200 

764.59683 

750.51859 

736.39471 

722.23292 

708.03501 

20.400 

693.80334 

679.53811 

665.23596 

650.90422 

636.54602 

20.600 

622.15676 

607.74401 

593.31063 

578.85846 

564.391-3 

20.800 

549.91232 

535.42714 

520.93674 

506.44745 

491.96713 

21.000 

477.49537 

463.03894 

448.60418 

434.20319 

419.83668 

21.200 

405.51304 

391.25183 

377.04965 

362.91724 

348.86909 

21.400 

334.92293 

321.07700 

307.35273 

293.76725 

280.32244 

21.600 

267.03443 

253.91885 

241.03571 

228.33303 

215.84614 

21.800 

203.61188 

191.69625 

180.01386 

168.58645 

157.48347 

22.000 

146.74482 

136.37370 

126.42682 

116.82968 

107.71723 

22.200 

98.99905 

90.77354 

82.98512 

75.56208 

68.60395 

22.400 

62.09665 

55.95399 

50.21239 

44.84255 

39.80813 

22.600 

35.14212 

30.80437 

26.82948 

23.09533 

19.62218 

22.800 

16.43807 

13.48717 
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Figure  B2  : Comparison  of  experimental  data  with  predated  spectra  using  assumed 
concentration  profile  1 for  sample  S-5-3  at  750  KeV. 


Figure  B6  Comparison  of  experimental  data  with  predicted  spectra  using  assumed 
concentration  profile  2 for  sample  S-5-3  at  800  KeV. 


Figure  B8  :Compar1son  of  experimental  data  with  predicted  spectra  using  assumed 
concentration  profile  3 for  sample  $-5-3  at  750  KeV. 


y 


Figure  BlO-.Assumed  concentration  profile 
4 for  sample  S-5-3. 
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Figure  B11  rComparlson  of  experimental  data  with  predicted  spectra  using  assumed 
concentration  profile  4 for  sample  S-5-3  at  750  KeV. 


Figure  B1 2: Comparison  of  experimental  data  with  predicted  spectra  using  assumed 
concentration  profile  4 for  sample  S-5-3  at  800  KeV. 


Figure  B1 4:  Comparison  of  experimental  data  with  predicted  spectra  uslngassumed 
concentration  profile  5 for  sample  S-5-3  at  750  KeV. 


Figure  B1 5 '.Comparison  of  experimental  data  with  predicted  spectra  using  assumed 
concentration  profile  5 for  sample  $-5-3  at  800  KeV. 
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Channel  Number 

Figure  : Fit  of  Sample  S-4-1  at  an  incident  proton  energy  of  750  KeV. 


Figure  C^:  Fit  of  Sample  S-5-1  at  an  Incident  proton  energy  of  750  KeV 
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